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I. INTRODUCTION 

Beams of atoms and molecules are stalwart tools for 
spectroscopy and studies of coUisional processes [ij [5]. 
The supersonic expansion technique can create cold 
beams of many species of atoms and molecules. How- 
ever, the resulting beam is typically moving at a speed 
of 300-600 m s~^ in the lab frame, and for a large class 
of species has insufficient flux (i.e. brightness) for im- 
portant applications. In contrast, buffer gas beams[31 H] 
can be a superior method in many cases, producing cold 
and relatively slow molecules in the lab frame with high 
brightness and great versatility. There are basic differ- 
ences between supersonic and buffer gas cooled beams 
regarding particular technological advantages and con- 
straints. At present, it is clear that not all of the possible 
variations on the buffer gas method have been studied. 
In this review, we will present a survey of the current 
state of the art in buffer gas beams, and explore some 
of the possible future directions that these new methods 
might take. 

Compared to supersonic expansion, the buffer gas 
cooled beam method has a fundamentally different ap- 
proach to cooling molecules into the kelvin regime. 
The production of cold molecules (starting from hot 
molecules) is achieved by initially mixing two gases in 
a cold cell (with dimensions of typically a few cm). The 
two gases are the hot "species of interest" molecules (in- 
troduced at an initial temperature Tq typically between 
300-10,000 K) and cold, inert "buffer" gas atoms (cooled 
to 2-20 K by the cold cell). The buffer gas in the cell is 
kept at a specifically tuned atom number density, typi- 
cally n — 10^^"^'' cm~^, which is low enough to prevent 
simple three body collision cluster formation involving 
the target molecule, yet high enough to provide enough 
collisions for thermalization before the molecules touch 
the walls of the cold cell. A beam of cold molecules 
can be formed when the buffer gas and target molecules 
escape the cell through a few-millimeter-sized orifice, 
or a more complicated exit structure, into a high vac- 
uum region as shown in Figure [2] . For certain buffer 
gas densities, the buffer gas aids in the extraction of 
molecules into the beam via a process called "hydrody- 
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namic enhancement." (S]. Finally, in the case where the 
mass of the molecule is higher than that of the buffer 
gas atom, there is a velocity-induced angular narrowing 
of the molecular beam, which increases the on-axis beam 
intensity ,5J. Although such enhancement has long been 
recognized in room-temperature beams [6], it is seldom 
employed because it requires intermediate Reynolds num- 
bers, in conflict with the high Reynolds numbers neces- 
sary for full supersonic cooling of molecules in the beam. 
In buffer gas cooled beams, on the other hand, the inter- 
mediate Reynolds number regime is often ideal for cool- 
ing, and allows this technique to take full advantage of 
the intensity enhancement from angular narrowing. 

With cryogenic cooling, high gas densities are not 
needed in the buffer cell to cool into the kelvin regime. 
In the case of supersonic beams, the high gas densities 
required in the source can be undesirable. In the case 
of buffer gas beams, the cryogenic environment and rel- 
atively low flow of buffer gas into the high vacuum beam 
region allows for 100% duty cycle (continuous) beam 
operation, without relying on external vacuum pumps. 
Rather, internal cryopumping provides excellent vacuum 
in the beam region. This combination of characteristics 
has led to high-brightness cold molecule sources (see Ta- 
bles |l] and |ll]) , for both chemically reactive (e.g. pulsed 
cold ThO, producing 3 x 10^^ ground state molecules 
per steradian during a few ms long pulse j7|) and stable 
molecules (e.g. continuous cold O2, producing w 3 x 10^'^ 
cold molecules per second [51 [H]). 



A. Cold atoms and molecules 

In the past two decades, the evaporative cooling of 
atoms to ultracold temperatures at high phase space 
density has opened new chapters for physics and led to 
exciting discoveries, including the realization of Bose- 
Einstein condensation|20|, |2T] . strongly correlated sys- 
tems in dilute gases [lH US], and controlled quantum 
simulation [25] . Meanwhile, the success of cold atom 
methods and new theory has inspired the vigorous pur- 
suit of molecule cooling. Molecules are more complicated 
than atoms and possess two key features not present in 
atoms: additional internal degrees of freedom, in the 
form of molecular rotation and vibration, and the pos- 
sibility (with polar molecules) to exhibit an atomic unit 
of electric dipole moment in the lab frame, which can 
lead to systems with long range, anisotropic, and tunable 
interactions. The rich internal structure and chemical di- 
versity of molecules could provide platforms for explor- 
ing science in diverse fields, ranging from fundamental 
physics, cold chemistry, molecular physics, and quantum 
physics[26]. We will list just a few of these applications 
here. 

• Molecules have enhanced sensitivity (as compared 
to atoms) to violations of fundamental symmetries, 
such as the possible existence of the electron electric 
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FIG. 1. Schematic velocity distributions for selected effusive, 
supersonic, and buffer gas cooled beam sources. The buffer 
gas cooled beam properties are taken from a ThO source with 
neon buffer gas [7]. The effusive beam is a simulated room 
temperature source of a species with mass of 100 amu, and 
the simulated supersonic source uses room temperature xenon 
as the carrier gas. Compared to the effusive beam, the buffer 
gas beam has a much lower temperature (i.e. smaller veloc- 
ity spread). Compared to the supersonic beam, the buffer 
gas beam has a comparable temperature but lower forward 
velocity. Supersonic sources typically have much higher aver- 
age forward velocity than the one presented above (~ 600 m 
s~^ for room temperature argon, or ~ 1800 m s~^ for room 
temperature helium), [1] and effusive sources for many species 
(like ThO O [14]) would require oven temperatures of > 1000 
K, making the distribution much wider with a much higher 
mean. The distributions above are normalized; however, the 
buffer gas cooled source of many species has a considerably 
higher flux. See Table|l]for data about experimentally realized 
buffer gas cooled beams. 

dipole moment |27l I28j , and parity- violating nuclear 
moments [m ISO]- 

• The internal degrees of freedom of polar molecules 
have been proposed as qubits for quantum 
computers [31], and are ideal for storage of quan- 
tum information [26j |32l [33] . 

• The long-range electric dipole-dipolc interaction 
between polar molecules may give rise to novel 
quantum svstems [34l I35j . 

• Precision spectroscopy performed on vibrational or 
hyperfine states of cold molecules can probe the 
time variation of fundamental constants, such as 
the electron-to-proton mass ratio and the fine struc- 
ture constant [55H55] . 

• Studies of cold molecular chemistry in the labora- 
tory play an important role in understanding gas- 
phase chemistry of interstellar clouds, which can be 
as cold as 10 K [MSI]- 
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Method Species Intensity [s~^] Velocity [m s~^] 



Chemically reactive polar molecules 



Buffer gas 


ThO[3 


3 X 10" sr~^ 


170 


Buffer gas 


SrF[9] 


1.7 X 10^^ sr"^ 


140 


Buffer gas 


CaH[TD| 


5 X 10'' sr"^ 


40 


Supersonic 


YbF[n] 


1.4 X 10^'' sr"^ 


290 


Supersonic 


BaFp] 


1.3 X 10^" sr"^ 


500 


Effusive 


SrF[Tl] 


5 X 10" sr"^ 


650 


Effusive 


ThOlZlITl] 


1 X 10" sr-i 


540 



Stable molecule with significant vapor pressure at 300 K 
Buffer gas NDgiHIie] 1 x lO" 65 

Supersonic ND3[I71[T8] 1 x 10* 280 

Stark Decelerated ND3[IZ1[TH] 1 x 10® 13 
Effusive ND3[I9] 2 x 10^° 40 

TABLE I. A comparison of supersonic, effusive, and buffer gas cooled beams for selected molecules. The intensity is the number 
of molecules per quantum state per second, and the velocity is reported in the forward direction. If the angular spread of the 
beam is known, then the intensity per unit solid angle (i.e. brightness) is given. The ND3 sources are velocity-selected. For 
pulsed sources, the intensity is averaged over a time period that includes multiple pulses. For chemically reactive molecules, 
buffer gas sources can provide dramatically higher brightness, and with a slower forward velocity. In addition to brightness and 
forward velocity, we shall see later in this article that other features of buffer gas beams also compare favorably to supersonic 
and effusive beams. See Table |ll] for more data about experimentally realized buffer gas cooled beams. 



• Ultracold chemical reactions have been observed at 
a temperature of a few hundred nK, with reaction 
rates controllable by external electric fields 143 j . 

• Molecular collisions in the few partial wave regime 
reveal the molecular interaction in great detail [445- 

The most common (but not the only [17]) way to cool 
atoms to ultracold temperatures, defined as the temper- 
ature (typically <1 mK) where only s— wave collisions 
occur (for bosons), is laser cooling |48j. This technique 
relies on continuously scattering photons from atoms to 
dissipate the atom's motional energy. Typically ~ 10'' 
photon scattering events are needed to significantly re- 
duce the kinetic energy of the atom. Molecules gen- 
erally lack closed transitions that can easily cycle this 
many photons because the excited states of the molecules 
can decay to many vibrational or rotational states. Be- 
cause laser cooling molecules is more difficult than it is 
with atoms (and has only been demonstrated relatively 
recently [491150]). there continue to be broad efforts in de- 
veloping new molecule cooling methods in order to fully 
control the internal and motional degrees of freedom of 
molecules [26l [511153] . Additionally, many of these new 
proposed methods could work well with laser cooling, in 
some form. 

In general, cooling techniques for molecules can be 
broadly classified into two types: indirect and direct [5T1 
[55] . Indirect cooling relies on assembling two laser-cooled 
atoms into a bound molecule using photoassociation or 
magnetoassociation. The resulting molecules have the 
same translational temperature as their ultracold parent 
atoms but are typically in a highly excited vibrational 



state, which can be transferred (typically by a coher- 
ent transfer method, such as STIRAP|54]) to the abso- 
lute ground state. While indirect methods can access 
ultracold polar molecules with a high phase space den- 
sity, these molecules are currently limited to a combina- 
tion of the small subset of atoms which have been laser 
and evaporatively cooled, such as the alkali and alkaline 
earth atoms. Currently, high-density samples of ultra- 
cold polar molecules have only been demonstrated with 
a single species, KRb[55j. Since many atoms are not eas- 
ily amenable to laser cooling, there is a large class of 
molecules that are the focus of direct cooling; this class 
includes many of the molecules that are desirable for the 
applications listed above [26]. In this paper we shall focus 
on beams of molecules that are cooled directly through 
collisions with cryogenically cooled atoms. 



B. Buffer gas cooling and beam production 

At the heart of the buffer gas beam technique is 
buffer gas cooling, a direct cooling technique used in the 
first magnetic trapping of polar molecules over a decade 
ago 56]. The buffer gas cooling technique works by dis- 
sipating the energy of the species of interest via elastic 
collisions with cold, inert gas atoms, such as helium or 
neon. Since this cooling mechanism does not depend on 
the internal structure of the species (unlike laser cooling), 
buffer gas cooling can be applied to nearly any atom or 
small molecule [4], and certain large molecules [57]. He- 
lium maintains a sufficient vapor pressure down to a few 
hundred mK[¥[[55]. and the typical helium- molecule elas- 
tic cross section]!] of ~ lO""''* cm^ allows buffer gas cool- 
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FIG. 2. A schematic of a buffer gas beam cell, which is main- 
tained at a temperature of few K. (a) The buffer gas (typically 
helium or neon) enters the cell through a fill line and exits via 
a cell aperture on the other side of the cell. The important 
physical dimensions (cell length Lceii, cell diameter dceii, and 
aperture diameter daperture) are indicated, (b) Introduction 
of species via laser ablation, (c) Molecules thermalize to the 
buffer gas and form a molecular beam with the fiowing buffer 
gas. 



ing and trap loading to be realized with modest cell sizes 
(few X few X few cm"^). Buffer gas cooling (using both 
helium and neon) has been used to create beams, load 
magnetic traps, and perform spectroscopy in a cold buffer 
gas celH4J. In this paper we shall focus only on one as- 
pect, the creation (and characterization) of buffer gas 
cooled beams. 

As shown in Figure [2j a buffer gas beam is formed 
simply by adding an exit aperture to one side of a cold 
cell filled with a buffer gas. A constant buffer gas den- 
sity is maintained by continuously flowing the buffer gas 
through a fill line on the other side of the cell. Figure 
^h) shows laser ablation, one of several methods used to 
introduce molecules of interest into the cell (see discus- 
sion in Section II A 2 1. After production and injection of 



molecules in the buffer gas cell, the molecules thermalize 
with the buffer gas translationally and rotationally, and 
are carried out of the cell with the flowing buffer gas, 
forming a molecular beam (Figure [2]jc)). 

Cold molecular beams created from a buffer gas cell 
are made in a few key steps: the introduction of "hot" 
molecules in a "cold" buffer gas cell, the thermalization 
of the molecules with the buffer gas, and the extraction 
of molecules into a beam. Each of these steps affects 
the properties of the resultant beam. The detailed char- 
acteristics (forward velocity, velocity spreads, tempera- 
tures, and fluxes) of buffer gas beams in several operating 
regimes will be described in the next section. Techniques 
to manipulate the beam properties, such as guiding, fil- 
tering, cooling by isentropic expansion, and advanced cell 
geometries will also be discussed. The last part of this 
paper will discuss the applications of buffer gas beams, 
including the search for the electron electric dipole mo- 
ment using ThOt59j, the study of cold dipolar coUisions 
with NDajBO], and the realization of direct laser cooling 
of molecules with SrF [igil5D] . 



II. BUFFER GAS COOLED BEAMS 

The details of buffer gas cooled beam production and 
properties will be presented in this section. In our treat- 
ment, we shall assume that the buffer gas is a noble gas, 
and that the species of interest is seeded in the buffer gas 
flow with a low fractional concentration, about 1% or 
less. We will often refer to the species of interest as "the 
molecule," even though buffer gas cooled beams of atoms 
(e.g. Yb) are also of interest. As is typically the case, we 
shall assume that the species of interest is heavier than 
the buffer gas (though the analysis is easily extended to 
lighter molecules). 



A. Species production, thermalization, diffusion, 
and extraction 

In this section we present estimates of physical param- 
eters that can be used to support an intuitive under- 
standing of the processes occurring in the buffer gas cell. 
These derivations will all be approximate, and will vary 
depending on geometry, species, introduction method, 
temperature range, density range, etc. However, they 
are typically descriptive within an order of unity. 



1. Buffer gas flow through the cell 

Consider a buffer gas cell as depicted in Figure [2j The 
cell has a volume of Vceii = Aceii x Lceii, where Lceii is the 
length of the cell interior, and Aceii « (ice« cross- 
sectional area (which may be round or square, but has a 
characteristic length of dceii)- L^eii is typically a few cm, 
and Aceii is typically a few cm^. The cell is held at a fixed 
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temperature Tq by a cryogenic refrigerator, typically be- 
tween 1 K and 20 K. Buffer gas of mass rrib is introduced 
into the cell volume by a long, thin tube, or "fill line". 
The buffer gas exits the cell through an aperture of char- 
acteristic length daperture and area Aaperture (for the case 



where the aperture is a rectangle, d, 



aperture 



is the shorter 



dimension). Typical values are daperture = 1 — 5 mm, 
and Aaperture = 5 — 25 mm^ . A buffer gas flow rate fo,b 
into the cell can be set with a mass flow controller. Here 
the subscript "6" refers to the buffer gas, and "0" refers 
to conditions in the cell at steady-state, or "stagnation" 
conditions. The most commonly used unit for flow is the 
standard cubic centimeter per minute, or SCCM, which 
equals approximately 4.5 x 10^^ gas atoms per second. 
Typical flow rates for buffer gas beams are /o_f, = 1 — 100 
SCCM. Technical details of cell construction, and con- 
struction of a buffer gas cooled beam apparatus in gen- 
eral, are discussed in Section p^I E[ 

At steady state, the flow rate out of the cell is given 
by the molecular conductance of the aperture, font = 
jno^bVo^b Aaperture, wherc 



Wo, 6 



iksTo 
nmb 



(1) 



is the mean thermal velocity of the buffer gas inside the 
cell (about 140 m s^^ for 4 K helium or 17 K neon), kg is 
Boltzmann's constant, and no,b is the stagnation number 
density of buffer gas atoms. Therefore, the number den- 
sity no.b is set by controlling the flow via the steady-state 
relationship font = /o,6, or 



no,b 



4/o.i 



A 



aperture'^0,b 



(2) 



With typical cell aperture sizes and temperatures, a flow 
of 1 SCCM corresponds to a stagnation density of about 
10^^ cm^'^, so typical values for the stagnation density are 
10^^ — 10^^ cm~'^. Note that in the above equation we as- 
sume that the flow through the aperture is purely molec- 
ular. At higher number densities the flow can become 
more fluid-like, and then the flow rate out will change; 
however, the difference is about a factor of two [HI], so 
the above equation is suitable for approximation. 



2. Introduction of species 

The species of interest can be introduced into the buffer 
gas cell through a number of methods [4] , including laser 
ablation, laser- induced atomic desorption (LIAD), beam 
injection, capillary filling, and discharge etching. The 
most commonly used techniques for creation of buffer 
gas cooled beams are laser ablation and capillary filling, 
so we will focus here on those two methods. 

Laser ablation. In laser ablation (see Figure [2]), a high 
energy pulsed laser is focused onto a solid precursor tar- 
get. After receiving the laser pulse, the solid precur- 
sor can eject gas-phase atoms or molecules of the de- 
sired species, often along with other detritus. The actual 



mechanism for how gas phase species results from the ab- 
lation of the solid precursor is not simple, and depends 
on the relationship between the length of the laser pulse, 
and the time constants for electronic and lattice heating 
[HU [SS]. While pulsed ablation has been studied with 
pulse widths from femtoseconds to milliseconds, the most 
common ablation laser (used for the majority of experi- 
ments in Table [Tl|) is a pulsed Nd:YAG, which can easily 
deliver up to 100 mJ of energy in a few ns. Either the 
fundamental (1064 nm) or first harmonic (532 nm) modes 
are used, and neither seems to have a distinct advantage 
over the other, except for certain technical conveniences 
afforded by a visible laser. 

Vapors of metals such as Na [3] or Yb [S] can be easily 
produced by ablation of the solid metal; however, cre- 
ation of gas phase molecules by laser ablation can be 
more complicated. Molecules with a stable solid phase 
(such as Pb0^3 ) can be created by simply using the solid 
phase as a precursor, but unstable molecules require care- 
ful choice of precursor. It is often the case that the de- 
sired diatomic molecule MX has a stable solid form M^Xf, 
which is a glass or ceramic, as is the case with BaF, CaH, 
SrF, YO, and ThO (and others). In such cases, experi- 
ence has shown that pressing (and sometimes sintering) 
a very fine powder of the stable solid often yields the best 
ablation target [71 [5] . 

The most important advantage of laser ablation is that 
it can be used to create a very wide variety of atoms 
and molecules with high flux [3]. The main drawback is 
that ablation is typically a "violent," non-thermal pro- 
cess that can result in unusual behavior of the gases in 
the buffer cell, and in the resulting beam. Ablation tends 
to give fairly inconsistent yields [71 [S], and can create plas- 
mas and complex plumes [M] with temperatures of several 
thousand Kf65j. However, buffer gas cells ca n be de signed 
to allow proper thermalization (see Section HAS) of the 



species, which can mitigate many of these problems. 

Capillary filling. If the species of interest has an ap- 
preciable vapor pressure at convenient temperatures, one 
can simply flow the species into the cell through a gas fill 
line, just like the buffer gas [66] ■ This method is simple 
in principle, but in practice it introduces some technical 
challenges. The cell and fill line must be carefully engi- 
neered to keep the fill line warm while keeping the cell 
cold and preventing freezing of the species at the point 
where the fill line enters the cell. For species with an 
appreciable vapor pressure at room temperature, such as 
O2 0, a properly-designed fill line can be attached di- 
rectly to the buffer gas cell to flow the species. For species 
that require higher temperatures to obtain an apprecia- 
ble vapor pressure, for example a 600 K oven to produce 
atomic K vapor, a complex, multi-stage cell must be con- 
structed [67] to keep heat from the oven out of the cell, 
as shown in Figure [3j 

The primary advantage of the capillary filling method 
is that it can be used to create high-flux beams that are 
continuous and robust. The main drawback is the limited 
number of species of interest that have appreciable vapor 
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FIG. 3. Left: simple capillary filling scheme for loading a molecule (large red circles) with vapor pressure at 300 K, such as 
O2 [5] or ND3 [151 167] . The species fiows down a warm gas line and is mixed with the buffer gas (small blue circles) , which flows 
down a cold gas line. Right: Scheme for capillary injection of K atoms (large red circles) from a 550 K oven into a 15 K neon 
(small blue circles) buffer gas cell. Unlike the simple capillary scheme, K atoms cannot simply flow down a warm gas line into 
the buffer gas cell, as the thermal conduction and blackbody heat loads would heat the cell too much. Instead, a multi-stage 
cell with thermal standoffs heat links at each stage must be used to reduce the heat load on the cell to acceptable levels. The 
amount of neon buffer gas introduced into the intermediate temperature regions must be sufficiently large to ensure that the 
K atoms are entrained in the flow and make it into the cold cell. Reproduced from [Patterson, D., Rasmussen, J., & Doyle, J. 
M., New J. Phys. 11(5), 055018 (2009) http://dx.doi.Org/10.1088/1367-2630/ll/5/055018, with permission from the Institute 
of Physics. 



pressure at temperatures low enough so that the method 
is feasible. As the required temperature increases, the 
technical challenges rapidly multiply (see Figure [s]). 

It should be noted that regardless of technique, the 
density of the species is typically (though not in the case 
of some capillary filling schemes [TSl [TC] ) < 1% of the 
number density of the buffer gas. This fact will be im- 
portant later on, because it allows us to treat the the gas 
flow properties as being determined solely by the buffer 
gas, with the species as a trace component. 



3. Thermalization 

Before the species flows out of the cell, it must undergo 
enough collisions with the cold buffer gas to become ther- 
malized to the cell temperature. A simple estimate of the 
necessary number of collisions can be obtained by ap- 
proximating the species and buffer gases as hard spheres 
[551 EH] . The mean loss in kinetic energy of the species 
per collision with a buffer gas atom results in a mean 
temperature change of 



tion as a differential equation: 



ATg — — (Tg — Tl,)/ K, where k, \ 



[mb + msf 



(3) 



Here T denotes temperature, m denotes mass, and the 
subscripts "6" and "s" refer to the buffer gas and species 
of interest, respectively. Therefore, the temperature of 
the species after M collisions, Ts{Af), will vary as 

T,{N) - TsiAf - 1) = -(TsiAf - 1) - n)/K (4) 

If we treat Af as large and the temperature change per 
collision as small, we can approximate this discrete equa- 



dAf 



= -{T,{M)-n)/K. 



(5) 



Solving this differential equation yields the ratio between 
the species and buffer gas temperatures: 



Tb 



= 1 



1 



Tb 
Ts{0) 

-tT' 



- 1 



-M/k. 



(6) 



(7) 



where in the last equality we have assumed that the 
species is introduced at a temperature much larger than 
the cell (and therefore buffer gas) temperature, i.e. 
Ts{Q) 3> Tb- If we estimate mb ^ 10 amu, to^ ^ 100 
amu, Tb - 10 K, and T,{Q) - 1000 K, the species should 
be within a few percent of the buffer gas temperature af- 
ter ~ 50 collisions. Under typical cell geometry and flow 
conditions, the thermalization time is a few milliseconds, 
as shown in Figure [4j 

While this simple model is good enough to obtain an 
estimate of the thermalization, it is not exact. In ana- 
lyzing thermalization of ablation-loaded YbF molecules 
in a helium buffer gas, Skoff et al. |7D] found that the 
model above did not fit the data; however, replacing the 
constant buffer gas temperature with an initially larger 
temperature which then exponentially decays to the cell 
temperature resulted in good fits. The physical moti- 
vation for this model is that the buffer gas is initially 
heated by the ablation, and then cools as the buffer gas 
re-thermalizes with the cell walls. Skoff et al. present 
detailed analysis of thermalization in a buffer gas, and 
the reader is referred to their work for details. 



7 




1 .5 2 2.5 
time after ablation (ms) 

FIG. 4. Thermalization of CaH in a helium buffer gas 
cell of temperature ~ 2 KJTD]. (a) The absorption sig- 
nal from a single ablation pulse, (b) Temperature of the 
CaH molecules vs. time after ablation, determined by 
Doppler spectroscopy. Reproduced from [Lu, H. et al., 
Phys. Chem. Chem. Phys. 13, 18986-18990 (2011) 
http:/ /dx.doi. org/10. 1039/clcp21206k by permission of the 
PCCP Owner Societies. 



Thus the cryogenic cell should be designed such that 
the species experiences at least 100 collisions or so before 
exiting the cell. To find the thermalization length, i.e. 
the typical distance that the species will travel before 
being thermalized, we need the mean free path of the 
species in the buffer gas, which is given approximately 



A 



-6,0 



aperture'^O.b 



a/1 + ms/rrib Afo^bO-b~s\/ms/mb 



(8) 



where (Tb-s is the thermally averaged elastic collision 
cross section (since the cross section typically varies with 
temperature), we assume ^ m^, and have used Eq. 
^ in the second step. For typical values of crb-s ~ 10"^"* 
cm'^, this mean free path is ~ 0.1 mm. Therefore, the 
thermalization length for the species in the buffer gas 
cell is typically no more than 100 x 0.1 mm = 1 cm. 

Note that the above discussion has pertained only to 
translational temperatures, yet buffer gas cooling is also 
effective at thermalization of internal states. Typical ro- 
tational relaxation cross sections for molecules with he- 
lium buffer gas are typically of order arot 10^(^^^^^) 
cm^, which means that around ab-s/crot ^ 10 — 100 
collisions are required to relax (or "quench") a rota- 
tional state Since this is a typical number of col- 
lisions required for motional thermalization, buffer gas 
cooling can be used to make samples of molecules which 
are both translationally and rotationally cold. Vibra- 
tional relaxation is less efficient, since the cross sections 
for vibrational quenching are typically several order of 
magnitude smaller than those for translational or rota- 
tional relaxation |4|. Several experiments [9l [56l [72] have 
seen the vibrational degree of freedom not in thermal 
equilibrium with the rotational or translational degrees 



of freedom. Further discussion of internal relaxation of 
molecules may be found elsewhere [U |4j |6T]. 

4 ■ Diffusion 

Once the species of interest is introduced into the cell 
and thermalized, we must consider the diffusion of the 
species in the buffer gas. Understanding the diffusion 
is crucial since the buffer gas cell is typically kept at a 
temperature where the species of interest has essentially 
no vapor pressure, and if it is allowed to diffuse to the 
cell walls before exiting the cell it will freeze and be lost. 
The diffusion constant for the species diffusing into the 
buffer gas is[7T] 



D = 



f 27rfcBTo 



16(no,s + nQ^b)cb-s V M 



1/2 



(9) 



where ^ = m,smb/ {ms + mb) is the reduced mass. It 
should be noted that different works use different forms 
for the diffusion coefficient (4j, though they agree to 
within factors of order unity and are all suitable for mak- 
ing estimates. Making the approximations rio.s ^ «o,6 
and nis 2> nib, we find 



D = 



2TTkBTo 
nib 



1/2 



32 TiQ^b^b-s 



(10) 



After a time t, a species molecule will have a mean- 
squared displacement of 



97r 



«0,6 



16 rio.fcCTb-s 



t 



(11) 



from its starting point. [75] Since the characteristic length 
of the cell interior is the cross-sectional length dceiii we 
can define the diffusion timescale Tdiff as {l\x'^)(Tdiff) = 



Aceiu or 



16 AceiiriQ^b^b-s 



(12) 



The diffusion time is typically 1-10 ms. 

Skoff et al. [70] performed a detailed theoretical anal- 
ysis and compared the results to measured absorption 
images in order to understand diffusion of YbF and Li in 
a helium buffer gas. The diffusion behavior was studied 
as a function of both helium density and cell temperature 
for both species. At low densities they verified the linear 
relationship between Tdiff and n^^b from Eq. (|12|), and a 
departure from linearity at high density. While this was 
not the first time that this relationship was observed [74]. 
their data and analysis allows for an informative possi- 
ble explanation: At low buffer gas densities, the ablation 
ejects the YbF molecules so that they are distributed all 
across the cell, and the diffusion model discussed above 
is valid. At high density, the YbF molecules are local- 
ized near the target and therefore diffuse only through 
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higher-order modes with shorter timescales. Previous 
models for this behavior included formation of dimers 
or clusters[73]; however, the model proposed by Skoff et 
al. jTD] seems reasonable considering that this behavior 
has been seen for several species, including both atoms 
and molecules [Ml [TQl [74H76] 



5. Extraction from the bujfer cell 

So far we have not considered the beam at all, having 
focused entirely on the in-cell dynamics. Once the species 
of interest is created in the gas phase and cooled in the 
buffer gas cell, it is necessary that the species flow out 
of the cell so that it can create a beam. As discussed 
above, extraction of the species from the cell must occur 
faster than the diffusion timescale Tdiff, so let's work out 
an estimate of the extraction, or "pumpout" time. The 
rate at which the buffer gas out of the cell is given by the 
molecular conductance of the cell aperture. 



Nb = -;NbVo,bAaperture/Vcelh 



(13) 



where Nb indicates the total number of buffer gas atoms 
in the cell, and Nb is the rate at which they are flowing 
out of the cell. The solution is an exponential decay with 
timescale Tpump^ the pumpout time, given by 



T, 



pump _ * 



(14) 



aperture 



The pumpout time is typically around 1-10 ms. Note 
that the pumpout time also sets the duration of the 
molecular pulse in the case of a beam with pulsed load- 
ing. If the buffer gas density in the cell is high enough 
that the species of interest follows the buffer gas flow, 
then this is a good estimate for the pumpout time for the 
species of interest as well. We now define [S] a dimension- 
less parameter to characterize the extraction behavior of 
the cell 



"/cell 



' pump 



4 nQ^b<^b-sAaperture ^ CTfc-s/o^fc 
LcellVa.b 



97r 



-icell 



, (15) 



where in the last approximation we dropped the order- 
unity prefactor as this is simply an estimate. This pa- 
rameter ^ceii characterizes the extraction behavior, and 
can be divided into two limits (see Figure [5|. 

For ^ceii ^ li the diffusion to the walls is faster than 
the extraction from the cell, so the majority of species 
molecules will stick to the cell walls and be lost. This 
"diffusion limit" [5] is characterized by low output flux 
of the molecular species, and is typically accompanied 
by a velocity distribution in the beam that is similar 
to that inside the cell [3]. In this limit, increasing the 
flow (thereby increasing ^ceii) has the effect of increasing 
the extraction efficiency fexti defined as the fraction of 
molecules created in the cell which escape into the beam. 



The precise dependence of fext on ^ceii is highly variable, 
and has been observed to be approximately linear, expo- 
nential, or cubic for the various species which have been 
examined[3[51[3in]. 

For 7ceH ^ 1, the molecules are mostly extracted 
from the cell before sticking to the walls, resulting in a 
beam of increased brightness. This limit, called "hydro- 
dynamic entrainment" or "hydrodynamic enhancement" 
[5] , is characterized by high output flux of the molecular 
species, and is typically accompanied by velocity distri- 
bution which can vary considerably from that present 
inside the cell. In this regime, the extraction efficiency 
plateaus and can be as high as >40%[51in], but is typi- 
cally observed to be around 10%. 

It should be noted that while the parameter ^ceii can 
be used to estimate the extraction efficiency quite well 
[3 [71 in] , there have been instances where this simple esti- 
mate breaks down. According to Eq. ( 15 1, the parameter 
7ceij has no explicit dependence on the aperture diame- 
ter; however, Hutzler et al. [7] found that while varying 
the cell aperture diameter in situ without varying other 
parameters, the extraction efficiency of a ThO beam in a 
neon buffer gas began to decrease with decreasing aper- 
ture size, as shown in Figure [5| These measurements 
suggest that the cell aperture diameter should not be 
too small (< 3 mm) in order to achieve good extraction. 

Regardless of the value for ^ceii, thermalization (see 
Section II A 3 ) must occur on a timescale faster than ei- 



ther Tpump or Tdiff to cool the species. Since neither Tdiff 
nor Tpump impose strict constraints on the cell geometry, 
it is possible to have both good thermalization and effi- 
cient extraction, as demonstrated by the large number of 
high flux, cold beams created with the buffer gas method 
(see Table |ll|. 



B. Properties of buffer gas cooled beams 

Now that we have reviewed the requirements for intro- 
duction of the species, thermalization, and extraction, we 
may focus on the properties of the resulting beam. 



1. Characterization of gas flow regimes 

As mentioned earlier, the species of interest typically 
constitutes < 1% of the number density of the buffer 
gas, so the flow properties are determined entirely by 
the buffer gas. In this section we will introduce the 
Reynolds number, which we will use to characterize gas 
flow. Note that the Reynolds (and Knudsen) number 
refers exclusively to the buffer gas, and not the species 
of interest (for our treatment). For a more thorough dis- 
cussion about gas flow, see texts about beams and gas 
dynamics [UleillHT]. 

The Reynolds number is defined as the ratio of inertial 
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Species 


Output /Brightness 




I'll [m s ^ 


Au|| [m s 


Molecules 










BaFp] 


1.6 X 10" sr-^ pulse"^ 




- 


- 




5 - 500 X 10* sr"^ pulse" 


1 


40-95 


65 


CHsFlU]'^ 


— 




45 


35 


CF3H 


— 




40 


35 




- 




- 


- 


NDainZ]'^ 


3 - 200 X 10* 




60-150 


25-100 


NDallSlIlS]^ 


1 - 10 X 10^" s"^ 




65 


50 


NDalSO]'^ 


1 X 10" 




100 


- 


02[5]« 


3 X 10^^ s"^ 




- 


- 


PbOig 


3 - 100 X 10* sr"^ pulse" 


1 


40-80 


30-40 


SrF(9] 


1 - 12 X 10^" sr"^ pulse" 


1 


125-200 


60-80 


SrOiTT] 


3 - 100 X 10^ sr"^ pulse" 


1 


65-180 


35-50 


ThO[Z] 


1 - 100 X 10^" pulse"^ 




120-200 


30-45 


YbFlZOlES] 


- 




- 


- 


YOI75I 


- 




160 


- 


Atoms 










K[62]° 


1 X 10^^ sr-i S-' 




130 


120 


Na[3] 


2 - 400 X 10* sr"^ pulse" 


1 


80-135 


60-120 


RblHni'^'^ 


3 X 10^" s-i 




190 


25-30 


Yb[5] 


5 X 10" sr-i pulse-i 




130 




Ybia^ 


5 X 10^" pulse"^ 




35 





TABLE II. A list of molecules which have been cooled in buffer gas beams, along with any measured properties. The output is 
either reported as molecules per state per pulse for pulsed beams, or per second for continous beams. If an angular spread was 
measured, the brightness is reported as molecules per state per unit steradian (sr) per pulse or per second. U|| is the forward 
velocity and Any is the full-width at half-maximum of the forward velocity distribution. Beams are ablation loaded and pulsed 
unless otherwise noted. A dash (-) means that the beam property was not reported in the indicated references. Notes: A = 



slowing cell (see Section IIC2|, B = continuous beam, velocity selected or guided by electromagnetic fields, C — flux reported 



is after beam collimation, D = continuous beam, loaded via capillary. 



to viscous forces in a fluid flowf 



Re 



F, 



intertial 



2j2 



pw d 



pwd 



(16) 



where p is the density, w is the flow velocity, p is 
the (dynamic) viscosity, and d is a characteristic length 
scale. In terms of kinetic quantities, we make use of 
the relationship between mean free path, density, and 
viscosity jSlj 



1 



p 



pXv, 



(17) 



where v is the mean thermal velocity and A is the mean 
free path, to express the Reynolds number as 



Re 



pw 



IpXv/d J/d"^ iKn' 



(18) 



where w/v Ma is the Mach number, and Kn = X/d 
is the Knudscn number. The Mach number for a gas of 
atomic weight m is defined as Ma — w/c, where 



(19) 



is the speed of sound in the gas, and 7 is the specific heat 
ratio. The value of 7 for a monoatomic gas (the relevant 
case for buffer gas beams) is 7 = 5/3, so in this case 



0.8w, 



(20) 



and we are justified in approximating Ma ss w/v. This 
gives us the important von Karman relation [HI [83] 



Ma : 



1 



-KnRe. 



(21) 



To relate these quantities to the case under consider- 
ation, we must mention two important facts about gas 
flow from an aperture into a vacuum [T]. First, the 
most relevant geometrical length scale which governs the 
flow behavior is the aperture diameter, since the prop- 
erties of the beam are set almost entirely by collisions 
that occur near the aperture. Therefore the characteris- 
tic length scale appearing in the formula for Re should 



be d 



^aperture • 



Second, near the aperture the gas 



atoms are traveling near their mean thermal velocity, so 
Ma « 1. Combining these facts with Eqs. ^ and ([s]) 
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FIG. 5. Extraction from a buffer gas cell. Top: Frac- 
tion of Yb atoms extracted into a helium cooled buffer gas 
beam vs. the parameter ^ceii (Eq. (15l). Reproduced 
from [Patterson, D. & Doyle, J. M., J. Chem. Phys. 126, 
154307 (2007) http://dx.doi.Org/10.1063/l.2717178 . Copy- 
right (2007), American Institute of Physics. Bottom: Frac- 
tion of ThO molecules extracted into a helium (x) and 
neon (o) buffer gas vs. the cell aperture size with a 
fixed flow rate. Reproduced from [Hutzler, N. et al., 
Phys. Chem. Chem. Phys. 13, 18976-18985 (2011) 
http://dx.doi.org/10.1039/clcp20901a (2011)] by permission 
of the PCCP Owner Societies. 



number, and number of collisions near the aperture are 
all linearly related. The possible types of flow can be 
roughly divided into three Reynolds number regimes, 
each of which will be discussed in the remainder of this 
section: 

• Effusive regime, Re < 1: In this regime there are 
typically no collisions near the aperture, so the 
beam properties are simply a sampling of the ther- 
mal distribution present in the cell. This regime is 
discussed further in Section Hi D II 

• Intermediate, or partially hydrodynamic regime, 
1 ^ < 100: Here there are enough collisions 
near the aperture to change the beam properties 
from those present in the cell, but not enough so 
that the flow is fluid-like. Buffer gas beams typi- 
cally operate in this regime; however, we will see 
examples of buffer gas beams in all three regimes. 

• Fully hydrodynamic, or "supersonic" regime, 100 < 
Re: In this regime the buffer gas begins to behave 
more like a fluid, and the beam properties become 
similar to those of a beam cooled via supersonic ex- 
pansion. This regime is discussed further in Section 
HID 21 



2. Forward velocity 
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allows us to write the relevant expression for Re in our 
situation, 



Re 



2d 



■aperture 



8V2fo., 



bf^b-b 



A 



h-6,0 



^aperture^O ,b 



(22) 



We shall parameterize the beam by the Reynolds num- 
ber Re of the buffer gas flow at the aperture. Ac- 
cording to Eq. (22), Re is approximately 2Kn^^ = 
'^daperture/ ^b-b, Or twicc the number of collisions within 
one aperture diameter of the aperture, i.e. "near" the 



aperture. From Eqs. (22) and ^ we can see that the 
in-cell buffer gas density, buffer gas flow rate, Reynolds 



FIG. 6. A schematic representation of beam forward veloc- 
ity vs. Reynolds number. In the effusive regime {Re < 1), 
the forward velocity is the thermal velocity of the (heavy) 
species. In the intermediate regime, collisions of the species 
with the buffer gas near the aperture accelerate the species; 
the velocity increase is linear with the Reynolds number until 
Re « 10, and then begins to asymptote to the final value. In 
the supersonic regime {Re > 100), the species has been fully 
accelerated to the forward velocity of the (light) buffer gas. 
Each of these regimes is discussed in detail below. 

If the beam is in the effusive regime, then there are 
typically no collisions near the aperture, and the forward 
velocity of the beam vu g is given by the forward velocity 
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of an effusive beam (Eq. (47)), where the appropriate 



thermal velocity is that of the species, i.e. 



Stt. 



(Re < 1) 



(23) 



As we will discuss, creating a purely effusive buffer gas 
beam without using advanced cell geometries (Section 
II C 2 1 can be challenging. 



In the intermediate regime, the molecules undergo col- 
lisions with the buffer gas atoms near (i.e. within one 
aperture diameter of) the cell aperture, whose average 
velocity is vo,b- This is larger than that of the (typically) 
heavier species Wq.s by a factor of -y/ nig /mi,. Since the 
collisions near the aperture are primarily in the forward 
direction, the species can be accelerated, or "boosted," 
to a forward velocity, V||,s, which is larger than the ther- 
mal velocity of the molecules (just as with supersonic 
beams [T]). 

If there are few collisions, we can estimate the rela- 
tionship between Re and v\\ s with a simple model from 
Maxwell et al.[3]. Near the aperture, the molecules un- 
dergo approximately Re/2 collisions. Each of these col- 
lisions gives the molecules a momentum kick in the for- 
ward direction of about « mf,wt,, so the net velocity boost 
is ~ Vb'mhRe/2mmoi- Since there are a small number of 
collisions, the forward velocity of the buffer gas is approx- 
imately given by the mean forward velocity of an effusive 
beam (Eq. (|47)), v\\^h ~ l-2wo,6, so for 1 < i?e < 10 (we 
shall justify tEe Re < 10 cutoff later on), 



0.6vo bRe 



nib 



(24) 



Therefore, the forward velocity should increase linearly 
with Re (and therefore with in-cell buffer gas density, or 
buffer gas flow) . Several papers [3 |T| H] have considered 
the behavior of s vs. i?e, and have seen this linear de- 
pendence at low Re. The more recent papers [71121 show 
data where the Reynolds number is large enough that a 
departure from the linear regime can be seen. Hutzler et 
al. [7] calculated the slope of the wy^s vs. Re relationship 
at low flow and found good agreement with the above 
model. 

This model necessarily breaks down as .j approaches 
~ Wo, 6 1 since the maximum possible forward velocity is 
1-4'i'o.b (from Eq. (54|), the forward velocity of a fully 



hydrodynamic buffer gas expansion. We therefore expect 
that the forward velocity should saturate to this value at 
large enough Re. To model the behavior outside of the 
linear regime considered above, Hutzler et al. [7] used the 
"sudden freeze" model [ET] , where it is assumed that the 
species molecules are in equilibrium with the buffer gas 
until the point along the beam where the buffer density is 
decreased enough such that there are no more collisions 
and the beam properties stop changing, or "freeze." The 
functional form for this model is (for Re > 10) 



,(i?e) « 1.4wo,fcVl -4i?e-4/5^ 



(25) 



and fits the data fairly well[7]. The estimate for where 
the linear to sudden-freeze model transition occurs is by 
finding at which flow there are collisions at a distance 
larger than one aperture diameter from the aperture, 
which happens for i?e > 10 [7]. 

Finally, for large enough Re, there should be enough 
collisions to fully boost the molecules to the forward ve- 
locity of the buffer gas. 



1.45, 



'o,b 



(Re > 100) 



(26) 



where the cutoff Re > 100 means that U|| ^(100) « 95% x 



W|[ according to Eq. (25), and from experiment [71 [3]. 
This limit corresponds to the supersonic flow regime. 
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FIG. 7. Forward velocity of extracted beam vs. buffer 
gas flow rate. Top: Forward velocity for Na and PbO 
in a helium cooled buffer gas beam. Reproduced from 
[MaxweU et al., Phys. Rev. Lett. Vol. 95, 173201 (2005) 
,Iittp:// dx.doi.org/10.1103/PhysRevLett. 95. 173201 Copy- 
right (2005) by the American Physical Society. Bottom: 
Forward velocity for ThO in a neon cooled buffer gas beam. 
Here the beam starts in the linear regime, and then fully sat- 
urates at the neon supersonic forward velocity. Reproduced 
from [Hutzler, N. et al., Phys. Chem. Chem. Phys. 13, 
18976-18985 (2011) http://dx.doi.org/10.1039/clcp20901a 
(2011)] by permission of the PCCP Owner Societies. 
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3. Velocity spreads 

Forward, or longitudinal velocity spread. In the effu- 
sive regime, the longitudinal velocity spread is the spread 
(FWHM) of the thermal ID Maxwell-Boltzmann distri- 
bution, 



/81n(2)fcf,ro 



wo,.A/i"ln(2) w 1.5vo,s- (27) 



Maxwell et al.fH] operated a buffer gas beam of PbO 
in this regime and found that the longitudinal velocity 
spread (as well as the transverse velocity spread, and ro- 
tational level distribution) were all in agreement with a 
thermal distribution at the cell temperature. 

If the Reynolds number of the flow is high enough, then 
the forward velocity spread can begin to decrease due to 
isentropic expansion of the buffer gas into the vacuum 
region. The translational temperature in the longitudinal 
(i.e. forward) direction can be decreased below the cell 



temperature, as will be discussed in Section II C 1 

Transverse velocity spread. Similarly, the transverse 
spread in the effusive regime is 



•0,s- 



(28) 



Also similar to the case with forward velocity, in the 
intermediate regime there will be collisions between the 
species and buffer gas near the aperture which can in- 
crease the transverse velocity spread. For low Re, we can 
estimate the transverse velocity spread near the aperture 
with a that is similar to the model for forward velocity. 



{Re/2)vo.b-,^ 



nib 

■^aperture "^mol 



(29) 



Hutzler et al. [7] used this model to describe the behavior 
of the transverse spread of a ThO beam in neon buffer 
gas and found it approximated the ratio A?;j^.s/i?e well. 

Experimentally, we are less concerned with the trans- 
verse velocity spread near the aperture than with the 
final velocity spread downstream, after collisions have 
ceased. This is more complicated to model since it con- 
flates the dynamics discussed above with the dynamics 
of the expansion. However, a similar argument shows 
that the transverse spread should increase linearly (above 
some Re) with increasing Re. This was examined in 
detail[71 [3], and a linear relationship was seen, with a 
ratio Av±^s/Rs that was ^ 2 times the ratio near the 
aperture [7]. 

Finally, as Re is further increased, the transverse ve- 
locity spread should saturate at the transverse spread of 
the buffer gas. Barry et al.fSl started to see this satura- 
tion behavior for SrF in a helium buffer gas start around 
Re w 100. 

We have seen that the shape of the Re vs. Av±^s re- 
lationship (see Figure [s]) is very similar to that which 
is plotted in Figure [6] However, as we shall see in the 



following section, the Reynold's numbers where the tran- 
sitions occur do not have to be the same as in the case 
of forward velocity. 
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FIG. 8. Behavior of transverse velocity spread vs. Reynolds 
number. Top: transverse velocity spread of SrF in a he- 
lium buffer gas cooled beam vs. Re at a fixed distance of 
20 mm (~ 7daperture) from the cell aperture j^. At low Re 
the transverse spread is in the linear regime, but it begins 
to saturate around Re ~ 50. Reproduced from [Barry et 
al., Phys. Chem. Chem. Phys. 13, 18936-18947 (2011) 
http://dx.doi.org/10.1039/clcp20335e by permission of the 
PCCP Owner Societies. Bottom: transverse velocity spread 
of ThO in a neon buffer gas[7] vs. Re at a fixed distance 
of 1 mm (< daperture) from the cell aperture. The transi- 
tion from flat (and approximately equal to the spread from 
a thermal distribution at To) to linearly increasing occurs 
around Re « 10. The data is well described by the model 
mentioned in the text. Reproduced from [Hutzler, N. et 
al., Phys. Chem. Chem. Phys. 13, 18976-18985 (2011) 
http://dx.doi.org/10.1039/clcp20901a] by permission of the 
PCCP Owner Societies. 

Typically, the transverse velocity spread is not a di- 
rectly important parameter. Actual experiments per- 
formed with atomic and molecular beams are generally 
performed at a distance from the cell aperture that is 
many times larger than daperture, and only sample a very 
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small solid angle of the output beam. This is both to 
reduce the gas load on the vacuum apparatus (the un- 
used portions of the beam can be differentially pumped 
to maintain a good vacuum), and to select atoms or 
molecules with nearly-identical directions of flight to re- 
duce Doppler broadening of transversely probed spec- 
troscopic transitions. However, the transverse velocity 
spread is important because it factors in to the diver- 
gence of the molecular beam, which will be discussed in 
the next section. 



4- Angular spread and divergence 

As mentioned earlier, the experimentally useful por- 
tion of an atomic or molecular beam is the portion that 
makes it through a typically small detection region at a 
large distance away. For this reason, an important char- 
acteristic of a beam is the value of the angular density 
distribution of the molecular beam. A parameter that 
can be used to characterize the angular spread of the 
beam is the full- width at half-maximum (FWHM) of the 
angular distribution, A^. This parameter is defined by 
equating N{A9) — ^Ntotai, where N{9) is the density at 
an angle of 9 from the aperture normal at a fixed distance 
from the aperture (see Section II D 1 1 . This parameter is 



We can then also define the solid angle spread Afi as 
the solid angle subtended by A9, i.e. 



very simple to calculate from transverse and longitudinal 
Doppler spectroscopic data; if a beam has a transverse 
velocity spread Av± and forward velocity W|| , then (see 
Figure ^ 



tan(A6l/2) 



Avj_/2 



(30) 



Ae = 2 arctan (^^^;^) (31) 




FIG. 9. An illustration showing the relationship between the 
forward velocity uy, transverse velocity spread Av±, and an- 
gular spread Ad from Eq. (1311), tan(A6l/2) = (Ai;j_/2)/u|| . 



A17 = 27r(l - cos(A6'/2)) 



(32) 



Note that this definition ignores the fact that at each 
point in the beam there is a velocity spread due to 
the non-zero temperature of the beam, so the measured 
transverse Doppler spread is a convolution of both beam 
translational temperature and the actual shape of the 
beam. However, if we measure this spread far enough 
away from the aperture that collisions have stopped, so 
that the aperture can be regarded as a point source from 
which the molecules are ballistically expanding, then in- 
ferring a spatial spread from the velocity spread is a valid 
approximation. 

For a buffer gas beam in the effusive regime, the angu- 
lar spread is given by Eq. (481 as A9 — 120°, and solid 



angle spread is given by Afl = tt. 

In the range 1 ^ -Re < 10, the forward velocity begins 
to increase linearly with Re yet the transverse velocity re- 
mains constant at Av±^s = l-5wo,s from Eq. (28). There- 
fore, Eq. (31) tells us that the divergence wiir begin to 
decrease. As the molecules begin to be boosted to the 
forward velocity of the buffer gas W|| ,j ss f || f, sa vo.t, the 
divergence should approach 



A9 = 2 arctan 




(33) 
(34) 
(35) 



where in the last approximation we assumed that ^ 
mil. The solid angle spread is then approximated by 

Ail — 27r(l — cos(2^ mi,/ms)) ~ nmii/'ms. (36) 
Notice that this can be much smaller than the corre- 



The gray shaded region indicates the spatial extent of the 
molecular beam, with darker indicating higher density. 



sponding spreads of tt for an effusive beam (Eq. (49 )), or 
1.4 for a supersonic beam (Eq. ([56])). Both Maxwell et 
al.[3] (with helium buffer gas) and Hutzler et al.[3 (with 
neon buffer gas) have operated beams in this regime, and 
have been able to see the linearly decreasing divergence. 

As Re is increased to the point where the transverse 
spread begins to increase, the divergence will stop de- 
creasing. Hutzler et al. [7] were able to see this transition 
for a ThO beam in neon buffer gas, shown in Figure |10| 
With helium buffer gas, both Hutzler et al.fZ] and Barry 
et al.^9j observed that the divergence was constant over 
a combined range of Reynolds numbers 1 < Re < 150. 
This was due to the fact that the increases in transverse 
and forward velocities canceled each other almost exactly. 
This indicates that while there is some proposed "univer- 
sal shape" for the relationships of Re vs. wy and vs. Av±, 
the Reynolds numbers where transitions in behavior oc- 
cur for the two relationships need not overlap, and there 
is not a similar universal shape for Re vs. Af2. 
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FIG. 10. Divergence of a ThO beam in both a heUum and 
neon cooled buffer gas beam[7]. For low Reynolds num- 
ber with neon buffer gas, we are in the regime of linearly 
increasing forward velocity but constant transverse veloc- 
ity spread (see Figures [t] and [sjl. Eventually the forward 
velocity stops increasing and the transverse velocity starts 
increasing, which results in the divergence ceasing to de- 
crease. With helium buffer gas, the divergence is essentially 
constant, similar to what was seen with SrF in a helium 
buffer gas cooled beam[5]. Reproduced from [Hutzler, N. 
et al., Phys. Chem. Chem. Phys. 13, 18976-18985 (2011) 
|http://dx.doi.org/10.1039/clcp20901a (2011)] by permission 
of the PCCP Owner Societies. 



Relationship of flow and extraction regimes 



By examining Eq. ( 22 ) 



for the Reynolds number, 
for the 



(15 



which governs the gas flow regime, and Eq 
extraction parameter jceii , which governs the species ex- 
traction regime, we can see that they are related by a 
factor which depends on the geometry. 



^cell ^ ^aperture 

-Re Lneii 



(37) 



This means that, at least in principle, it should be pos- 
sible to design a buffer gas source which is any combi- 
nation of diffusive (low efficiency) or hydrodynamically 
extracted (high efficiency), and with effusive, intermedi- 
ate, or hydrodynamic flow. Most buffer gas sources op- 
erate either in the effusive or intermediate regimes, and 
it is experimentally challenging to design a beam which 
is completely effusive, has good extraction, and sufficient 
thermalization. A purely effusive beam should have a 
forward velocity which, according to Eq. (47), does not 



change with source pressure; however, in buffer gas beam 
sources with good extraction [3 [71 IHl [H] the forward ve- 
locity of the molecules indeed changes with source pres- 
sure. Slowing cells (Section |ll C 2 1 can offer near-effusive 
velocity distributions, but typically have ~ 1% extraction 
effjciencvplITU]. 



6. Choice of buffer gas 

For many applications, helium is a natural choice of 
buffer gas. Helium has large vapor pressure at 4 K (that 
corresponds to an atom number density of > 10^^ cm""^), 
which is a convenient temperature for cryogenics: Many 
refrigerators can cool to ^ 4 K, as can a simple liquid he- 
lium cryostat. Helium also has the advantage that it has 
large vapor pressure at even lower temperatures (< 1 K), 
where operating the cell will result in a colder and slower 
beam[5l I67j . Simple cryogenic techniques can achieve 
temperatures of ^ 2 K and still handle the necessary 
heat loads; lower temperatures require a more compli- 
cated setup. 

Neon has also been used as a buffer gas [H [67], and it 
has some advantages of its own. The thermal velocity of 
the buffer gas scales as ^jT\)lm\i (Eq. (Il|), so as long as 
the cell temperature is kept below 20 Kwith neon buffer 
gas (whose mass is ~ 20 amu), the ther mal ve locity (and 
therefore forward velocity; see Section nB2) should be 



comparable to that of a 4 K helium-4 buffer gas beam 
source. Since neon has large enough vapor pressure down 
to ^ 14 K (where the density drops below ^ 10^^ cm^'^), 
it can compare favorably to a 4 K helium source in terms 
of forward velocity. 

While the forward velocity of a neon or helium buffer 
gas cooled beam can be comparable, in certain situations 
the temperatures in the beam can be comparable as well. 
Hutzler et al. [7| found that operating a buffer gas cooled 
beam resulted in isentropic cooling (see Section 11 C 1 ) 
that reduced the temperature of both a helium and a 
neon cooled beam of ThO to similar temperatures (see 
Table |ll| , even though the cell temperature was around 
17 K for the neon cooled beam and around 4 K for the 
helium cooled beam. This means that the beam proper- 
ties of a helium or neon cooled beam can be similar, with 
the neon beam offering some distinct technical advan- 
tages. Neon can be efficiently cryopumped by a 4 K sur- 
face, while helium requires a large-area adsorbent such as 
activated charcoal, as discussed in Section [ll E 3| Char- 
coal cryopumps are very effective; however, they must 
be emptied periodically, have reduced pumping speed as 
they begin to fill, and can influence the resulting molec- 
ular beam properties [71 [S]. Cryopumping of neon onto 
a cold surface does not display these properties j7]. Ad- 
ditionally, the properties of pulsed beams with helium 
buffer gas can show time-dependence not present with 
neon mil]. 



C. Additional cooling and slowing 

As discussed in the preceding sections, the properties of 
a buffer gas beam can be tuned considerably by altering 
the flow and aperture sizes to control the parameters Re 
and 7ce/i- However some applications require even finer 
control; for example, collision studies may require a very 
narrow and tunable velocity distribution, and trap load- 
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FIG. 11. Plots comparing the forward velocity distributions 
using different buffer gasses|B7]. Top: Calculated distribu- 
tions for noble gases at a temperature where there is suffi- 
cient vapor pressure for beam operation. While helium has 
sufficient vapor pressure below 4.2 K (and even 1 K), 4.2 
K is a typical temperature for both cryogenic refrigerators 
and liquid helium cooled cryostats. While the neon distri- 
bution is broader (i.e. hotter), for applications where only 
the slowest molecules are used (such as in trap loading, see 
Section IIICl, neon has an advantage. Bottom: measured 
velocity distributions of beam of ND3 in helium and neon 
buffer gas at 12 K. If the cell must be operated at a higher 
temperature for technical reasons (for example, if there is 
a large heat load on the buffer cell from a capillary load- 
ing scheme), this data shows that neon has the advantage 
in terms of creating slow molecules. The data here is veloc- 
ity selected. Reproduced from [Patterson, D., Rasmussen, 
J., & Doyle, J. M., New J. Phys. 11(5), 055018 (2009) 
[http://dx.doi.org/10.1088/1367-2630/l l/5/0550l8| with per- 
mission from the Institute of Physics. 



ing typically requires a very slow distribution (see Section 
IIIC). In this section we shall discuss some methods to 
manipulate the output of a buffer gas beam, to either 



reduce temperatures, reduce forward velocities, or select 
velocity classes. These techniques help make buffer gas 
beams highly versatile, and allow the many applications 
discussed in Section IllTl 



The boosting effect discussed in Section II B 2 can be 
detrimental if one aims to produce a slow beam, which 
is one of the benefits of buffer gas beams over super- 
sonic beams and effusive "oven" beams. One would 
therefore prefer to operate the beam in the effusive flow 
regime, i.e. at a low Reynolds number (Eq. (22)). 



It is possible to keep the Reynolds number fairly low 
while maintaining good extraction {'yceii > 1) by using 
a slit-shaped aperture, as we will now discuss. Consider 
the aperture as having a short and long dimension, so 
Aaperture = dghort X diong- FoT & fixcd internal cell ge- 
ometry, buffer gas, species, and cell temperature, we can 
see from Eqs. (22 1 and (151 that 



Re cx dshortno,b 



~Re 



cx di 



ong 



(38) 
(39) 

(40) 



Therefore, by increasing diong but keeping Aaperture 
fixed, we can decrease Re without changing ^eeii or the 
buffer gas density. Note that simply changing the aper- 
ture size while leaving all other parameters fixed also 
has the effect of changing Re but not ^ceii (Eq. (15)); 
however, this is not ideal for two reasons: First, the in- 
cell buffer gas density is constrained to be large enough 
that the species is thermalized, but small enough that the 
molecules can diffuse away from the injection point [70]. 
and this method would change the buffer gas density. 
Second, as mentioned in Section [lI A 5[ reducing the aper- 
ture size can in fact have an effect on the extraction for 
small enough aperture sizes. For these reasons, earlier 
buffer gas beam papers tended to use a slit [21 aperture, 
typically around 1x5 mm. 

While performing buffer gas beam studies that re- 
quired high flow rates to sufficiently thermalize atoms 
from a hot oven, [57| [5^ it was quickly realized that the 
downside of high Re flows could also come with a ben- 
efit: isentropic cooling from the free expansion of a gas, 
similar to what happens in supersonic beams [I]. This ef- 
fect was first seen (though not published) in K and Rb 
atoms cooled by a neon buffer gas|H71 HU], and it was 
subsequently fully characterized with ThO in neon and 
helium buffer gas beams at cell temperatures of about 
17 K and 4 K (respectively) [7] , and simultaneously with 
SrF in helium buffer gas at a cell temperature of 3 K[9 . 
In each case, the molecules were found to cool rotation- 
ally and translationally below the temperature of the cell, 
and the beam properties had very good qualitative agree- 
ment. Selected properties of these beams are listed in 
Table n 

The cell apertures for these studies were either round 
or square, and the studies were performed between Re ~ 
1 — 150. In each case, the Reynolds numbers were high 
enough that the rotational temperatures appeared to sat- 
urate around 2 K and 1 K for ThO and SrF respectively, 
as shown in Figure 12 The translational temperatures 
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(in the forward direction) of ThO in neon buffer gas was 
also decreased befow the temperature of the cell, and is 
shown in Figure [12] However, the translational temper- 
atures for both ThO and SrF in helium buffer gas were 
complicated by time-dependence of the temperature af- 
ter the ablation pulse, though in each case it was reduced 
below the cell temperature. 

It is worth noting that even in these fully-boosted, high 
Reynolds number regimes, where the forward velocity is 
fully saturated (see Section II B 2 1 , the velocity still com- 
pares very favorably to supersonic or effusive beams (see 
Tables [l| and [TT|) . In the next section, we shall see that 
there are additional techniques to further reduce the for- 
ward velocity of a buffer gas cooled beam. 
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2. Slowing cells 

As mentioned in the previous section, the beam prop- 
erties can be manipulated by choice of aperture geometry. 
This control can be extended by using an advanced aper- 
ture geometry known as a "slowing cell" [51 [TOj , shown 
in Figure |13[ which is used to create a region of inter- 
mediate pressure between the cell and vacuum (beam) 
regions, and acts as an effusive source that can still uti- 
lize the high extraction of the main cell. 

The main, or "first," cell is created like a typical buffer 
gas cell, and has the same typical buffer gas flow rates, di- 
mensions, etc. A second, or slowing, cell is then attached 
to the aperture of the main cell. This second cell has a 
very large area through which the buffer gas can flow, 
so the steady-state pressure inside it is low, about 10% 
of the pressure in the main cell. This makes the typical 
mean free path on the order of a few mm, so the species 
will scatter only a few times while in the slowing cell. 
If the slowing cell is also kept cold, these collisions will 
not heat the molecules, but will have the effect of reduc- 
ing their boosted forward velocity closer to the thermal 
velocity (see Section II B 2 ) . 



Species 



To [K] v^i 



N [puis 



Yb, Boosted!!] 2.6 130 

Yb, Slowed[5] 2.6 35 

CaH, Boosted[Tn] 1.8 110 

CaH, SlowedUO] 1.8 65 

CaH, SlowedllO] 1.8 40 



5 X 10^^ 
5 X 10^" 
5 X 10^ 
1 X 10^ 
5 X 10* 



TABLE III. A comparison of boosted bufler gas beam sources 
(i.e. with no slowing cell), to those with a slowing cell. 
"Boosted" means that collisions near the aperture accelerate 
the beam above the thermal velocity, as discussed in Section 
|II B 2[ "Slowed" means a slowing cell was added to the main 
cell. 
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FIG. 12. Additional cooling from isentropic expansion in 
a buffer gas cooled ThO beamljl. The cell temperature is 
around 17 K for neon buffer gas, and around 4 K for helium 
buffer gas. Top: Longitudinal velocity spread (temperature) 
in a neon cooled buffer gas beam, with two different aperture 
sizes. Bottom: Final rotational temperature in the beam, for 
comparing two different aperture sizes with neon buffer gas, 
and a single aperture size for helium buffer gas. The rotational 
temperature with helium buffer gas showed essentially no de- 
pendence on the aperture size. Figures are reproduced from 
[Hutzler, N. e t al., Phys. Chem. C hem. Phys. 13, 18976- 
18985 (2011) |http://dx doi.org/10 .1039/clcp20901a (2011)] 
by permission of the PCCP Owner Societies. 



The PbO beam of Maxwell et al. had a forward velocity 
of ~ 40 m at low flows, comparable to those achieved 
with Yb and CaH using slowing cells, but with an extrac- 
tion fraction of fext < 10~^, compared to fe^t ~ 10~^ for 
the slowing cells. 



Notice that operating a beam in the effusive regime can 
also result in low forward velocities (see Section II B 2 1 ; 
however, as mentioned earlier, it is difficult to operate 
a fully effusive beam that also has significant extraction. 



3. Guiding and velocity filtering 

Although they may not technically constitute cooling 
or slowing, electric or magnetic guides or velocity filters 
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FIG. 13. A slowing cell, adapted from Lu et al.[Tn]. Here the 
species is represented by red dots, and the helium buffer gas is 
represented by blue shading, where darker shading indicates 
higher density. The density in the slowing cell is kept low by 
gaps and vents, so that the species suffers only a few collisions 
in the slowing cell. These collisions can reduce the forward 
velocity without heating, while sacrificing flux. The method 
for introducing the species is not indicated in this figure, but 
has been accomplished with both ablation[5l[lO] and capillary 
filling [5]. 



can be used to change the output of a buffer gas cooled 
beam to suit a particular application's, 15, 16, 67, 83] ■ 
Guiding the species of interest can be useful for a number 
of reasons. First, if the presence of the buffer gas will in- 
terfere with the experiment, for example while perform- 
ing collision studies, a guide can be used to transport 
the species far from the buffer gas source where there 
is no appreciable buffer gas pressurelHO]. Second, if the 
beam is to be studied in a room temperature appara- 
tus, for example to perform laser cooling 50 or pre- 
cision measurements [55] ■ the beam may have to travel 
a long distance to exit the cryogenic apparatus before 
passing into the experimental region. Without a guide, 
the molecule density will decrease with the square of the 
distance from the cell, which could result in significant 
reduction. An electrostatic or magnetic guide or lens can 
help reduce this loss from beam divergence. 

Magnetic B or electric £ fields can be used to guide 
or focus a molecular beam by creating a restoring force 
perpendicular to the motion of the molecule [1]. If the 
molecular state has a magnetic moment /xg, then the in- 
teraction (of a low-field seeking state) with a magnetic 
field leads to a potential energy U = If we assume 

that the molecule is traveling in a cylindrically symmet- 
ric field B(r\ where r is the distance from the axis of 
symmetry, then the field from a rt— pole configuration 
is \B[r')\ = ;Bor'^"'~^^/^, and the potential energy of the 
molecule is C/ = ^eSor^""^)/^. In addition to guiding, 
multipole fields can also be used to focus molecules by 
creating a linear restoring force {U oc r^), which we can 
see requires a hexapole (n — 6) field configuration. 

The situation is similar for polar molecules. If a 
molecular state has a permanent electric dipole mo- 



ment \i£ (as is the case for states with A > 0, such 
as n or A states[85]), then the potential energy of the 
(low-field seeking state) molecule in an n— pole field is 
U = /ifrfo'^^""^^''^- However, unlike the case with mag- 
netic dipole moments, polar molecules can be in states 
(such as E states [85]) which have an induced dipole mo- 
ments, leading to a quadratic shift in fields. In this case, 
the dipole moment (for small fields) is proportional to the 
electric field /.tf oc |£|, so potential energy of the molecule 
is oc r"~^. If we wish to focus polar molecules 

with a linear restoring force, we can see that a hexapole 
is required for states with a permanent dipole moment, 
or a quadrupole (n — 4) for states with an induced dipole 
moment. 

If the maximum electric field in the guide (typically 
~ 10 — 40 kV cm~^) is £max and the molecules have 
Stark shift W{E) in a field then the transverse depth 
of the guide is defined as W^ax — W{£max), which is 
typically ~ 1 Kxfcs. This sets the maximum transverse 
velocity spread in the guide v±^max to be 

1 



^2W„,ax/ms, (41) 
which is around 10-30 m s~^. 
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FIG. 14. Electrostatic guide for NDsd?]. The cell is kept 
at 15 K and contains neon buffer gas (small blue circles) and 
ND3 (large red circles) introduced from a capillary. As the 
molecules exit the cell, an "electrostatic gate" can be turned 
on to direct molecules into the guide. This guide will only 
admit slow molecules around the bend; faster molecules will 
simply overcome the potential barrier and escape the guide. 
Adapted from Patterson et al. [67) 

In addition to guiding or focusing to increase useful 
flux, electric or magnetic fields can be used to select par- 
ticular velocity ranges from the output of the molecular 
beam, which is useful for collision studies and trap load- 



ing experiments (see Section IIICl, as will be discussed. 



Molecules of a particular forward velocity range can be 
selected by rotating mechanical objects [SB]; however, us- 
ing fields is much less technically challenging, especially 
in a cryogenic environment. The first realization of this 
method 87J (and later realizations with buffer gas cooled 
beams[5t il^[T51[57] ) was to simply bend one of the guides 
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described above, similar to what is shown in Figure 14 
Molecules with high enough kinetic energies in the for- 
ward direction will simply pass over the potential bar- 
rier, and only the more slow-moving molecules will re- 
main in the guide. More specifically, for a guide of depth 
Wmax, guide radius p (i.e. distance from the center of the 
guide to the location of the maximum electric field), and 
guide radius of curvature R, the cutoff for longitudinal 
velocities which will remain in the guide is found[51] by 
equating the centrifugal force msfjj ^^^/^R to the restor- 
ing force ~ Wmax/p caused by the molecule traveling up 
the potential hill of length p, or 

^svl,„ax/R ~ Wmax/p ^ Wll 



\/WmaxRI prris 



(42) 

This technique can be used to filter only slow molecules 
from a pulsed or continuous buffer gas source, for exam- 
ple to load a trap. In fact, this technique can also be 
applied to thermal effusive sources; Rangwala et al.[87] 
filtered slow H2CO and ND3 molecules from a 300 K ef- 
fusive source and obtained distributions that had widths 
and means similar to those of a 5 K thermal distribution, 
while maintaining a high fiux of ss 10^ s^^ guided H2CO 
molecules. 

In addition to simply selecting slow molecules, elec- 
trostatic guides can be used to select a window of ve- 
locities by using multiple switched guides. Sommer 
et al.|84j devised a three-stage, switchable electrostatic 
guide, fed by a helium buffer gas cooled source of ND3 
molecules, which could deliver controllable, velocity- 
selected molecule pulses. By switching the individual 
guide segments on or off at precise times, both a lower 
and upper velocity cutoff can be imposed on molecules 
that remain in the guide. Velocity windows of FWHM 
^ 3 — 20 m s^^ wide, with centers between 20 and 100 
m s^^, could be efficiently selected with low losses, as 
shown in Figure [TS} 



D. Effusive and supersonic beam properties 

For the sake of comparison, we will briefiy review some 
properties of effusive and supersonic beams. Detailed 
discussions may be found in existing literature pi 161]. 
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FIG. 15. Velocity-selected pulses of ND3 obtained using mul- 
tiple bent, switchable guides 84'. The pulse centered about 
60 m s~^ contains approximately 10^ molecules. Reproduced 
from [Sommer, C. et al., Phys. Rev. A 82(1), 013410 (2010) 
|http://dx.do i.org/10.1 103 /PhysRev A.82.013410, . Copyright 
(2010) by the American Physical Society. 



of the species of interest is large at the source temper- 
ature (as opposed to buffer gas cooled effusive sources, 
which can operate at temperatures where the species has 
no appreciable vapor pressure). 

The number density in the beam resulting from a dif- 
ferential aperture area dA is given by 



n{R, V , 9) = n{R, 9) f {v) , where 
n{R,9) = ^^^^dA, 



f(v) = 



47ri?2 
32 



^2g-4«V7r(S)^ 



(43) 
(44) 

(45) 



R is the distance from the aperture, 9 is the angle from 
the aperture normal, hq is the stagnation density in the 
cell, n{R, 9) is the total number density distribution inte- 
grated over velocity, and f{v) is the normalized velocity 
distribution in the cell. The velocity distribution in the 
beam is given by |61| 



fbeamiv) = {v/v)f{v) 



32 



(46) 



1. Effusive beams 

In an effusive gas fiow from an aperture, the typical 
escaping gas atom has no collisions near the aperture 
{Kn > 1). Therefore, the resulting beam can be consid- 
ered as a sampling of the velocity distribution in the cell. 
A typical setup is a gas cell with a thin exit aperture, 
where the thickness of the aperture and the aperture di- 
ameter daperture are both much smaller than the mean 
free path of the gas at stagnation conditions. Note that 
the effusive beams under discussion in this section are 
oven-type effusive sources, i.e. where the vapor pressure 



From the velocity distribution we can extract the mean 
forward velocity of the beam. 



II, e// 



V fbeamiv) dv 
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V w 1.2ti 



(47) 



From the total number density in the beam n(R,9), 
we can extract the FWHM (full- width at half-maximum) 
of the characteristic angular spread A9 by solving 
n{R,M/2) = \n{R,Q), which gives 



M= — = 120°, 
3 



(48) 
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or the characteristic sohd angle Afi 

An = 27r(l - cos(A6'/2)) = tt 



(49) 



Note that this is the angular spread from a differential 
area of the aperture; near the aperture we would have 
to integrate over the area of the aperture to obtain the 
full shape, however this angular spread is valid in the far 
field. 

The rate at which molecules escape the cell and there- 
fore enter the beam is simply the molecular flow rate 
through an aperture. 



1 



N — -no.bVoAaperture- 



(50) 



Effusive beams of certain species have very large fluxes, 
such as certain metal atoms or low-reactivity molecules 
with high vapor pressure. For alkali metals [48], the cell 
(or oven) temperature required to achieve a vapor pres- 
sure of 1 torr is between around 500 and 1000 K. The flow 



through a 1 mm aperture is then around 10 
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Note that in this situation there is no buffer gas. 

The beam resulting from an effusive source is not im- 
mediately useful for many applications. The large for- 
ward velocity (typically several hundred m s~^) would 
limit laser interrogation time, and broad velocity dis- 
tributions would lead to significantly broadened spec- 
tral lines. Some atoms can be slowed and cooled using 
powerful optical techniques 148]: however, until recently 
molecules have (with few exceptions [IHl HO]) resisted 
these optical techniques due to their complex internal 
structure. In addition, molecules have internal degrees 
of freedom which can be excited by the large temper- 
atures in an oven. The rotational energy constant for 
diatomic molecules is typically 1-10 K xks, so at typical 
oven temperatures the molecules can be distributed over 
hundreds of rotational states. 



2. Fully hydrodynamic, or "supersonic" beams 

In a fully hydrodynamic, or "supersonic" beam, the 
gas experiences many collisions near the exit aperture 
Kn ^ 1 (typically Kn < 10~^), and therefore the beam 
properties are determined by the flow properties of the 
gas. In this case we cannot apply simple gas kinetics as 
in the case of effusive beams, but instead must consider 
the dymanics of a compressible fluid. We shall assume 
that the gas is monoatomic. 

A typical supersonic source has Pq ^ 1 atm, and 
daperture ~ 1 mm^ . Supcrsonic sources are often cooled 
so that the forward velocity is reduced, and typically 
have To — 200 — 300 K. Increasing the backing pres- 
sure Pq can be used to further reduce the temperature 
of the molecules, though at high enough pressures clus- 
ter formation can begin to negatively affect the beam 
properties [1 1 . On the other hand, buffer gas beams can 
be operated with high enough backing pressure such 



that increasing the pressure further will no longer reduce 
the beam temperature, without degradation of the beam 
properties [3 [5]. Additionally, for supersonic beams, the 
introduction of chemically reactive or refractory species is 
challenging due to the short mean free path between col- 
lisions in the source, and are typically introduced into the 
expansion plume pTl [55H90 j . This reduces the brightness 
of the beam. On the other hand, beams of polar radicals 
can be created with buffer gas cooling [3 IHj with around 
100 times the brightness of supersonic beams of similar 
species (see Table |l]). 

The gas flow rate from the aperture in a supersonic 
source can be on the order of 1 standard liter per second, 
which would make keeping good vacuum in the appara- 
tus difficult if the beam were to operate continuously. For 
this reason, supersonic beams are often pulsed to reduce 
time-averaged gas load on the vacuum system. Contin- 
uous, or "Campargue"-tvpeT1 WT\ supersonic beams are 
possible, although they introduce many technical chal- 
lenges. On the other hand, buffer gas cooled beams can 
be operated continuously without considerable difficulty 
(see Table |ll] for a list of continuous buffer gas cooled 
beams). 

For a monoatomic gas with specific heat ratio 7 — 
5/3, the number density and temperature in a supersonic 
expansion are related bv|61j 



n 



3/2 



(51) 



In the far field (typically more than four times the 
aperture diameter away from the aperture [HIDj the num- 
ber density will fall off as a point source, n{R) oc R~^, 
where R is the distance to the aperture. Therefore, 



n{R) oc R-^, T{R) oc R-^/^. 



(52) 



Thus we see that, unlike in the case of an effusive 
beam, as the beam expands the temperature decreases. 
The temperature will continue to decrease until the gas 
density becomes low enough that collisions stop, the gas 
ceases to act like a fiuid, and the gas atoms simply fly 
ballistically. This transition is often called "freezing" or 
"quitting." Even from a room-temperature supersonic 
source, it is not uncommon to have beam-frame temper- 
atures of around 1 K. 

The relationship between the forward velocity and 
temperature of an ideal monoatomic gas expansion is 
given by [61] 



(53) 



If this gas is allowed to expand a long enough distance 
such that T <C To, then the final forward velocity is 



SfcsTo 



Wo 



lAvQ. 



(54) 



20 



A standard supersonic source is argon expanding from 
a 300 K cell, which has a forward velocity of about 600 
m s^^. A more technically challenging source has argon 
expanding from a 210 K cell, which has a forward velocity 
of about 300 m s^^. On the other hand, buffer gas beams 
typically have a forward velocity of < 200 m s~^, a nd can 
be as slow as 40 m s~^ (see section Section IIB2|. 

The label "supersonic" comes from the fact that as the 
gas expands from the aperture, the temperature drops so 
that the speed of sound (Eq. (19l) decreases yet the for- 
ward (flow) velocity increases (Eq. (53)), so that even- 
tually the Mach number becomes larger than 1. In fact, 
setting the stagnation pressure and temperature so that 
the Mach number is exactly 1 at the aperture yields max- 
imum gas flow rate through the aperture 'BT . 

The number density as a function of the distance R 
and the angle 9 from the aperture normal is given ap- 
proximately by |92j 



n{R,e) = n{R,0)cos^ 



(55) 



where (/) w 1.4 for a monoatomic gas. We can then find 
the angular spread A9 and solid angle Afi as in Eqs. ( 48 ) 
and (|49| to be 



AO- 



1.4 ==79°, An^lA 



(56) 



If there is a small amount of a species of interest mixed 
in with the main (or carrier) gas, there are enough col- 
lisions that the species will follow the carrier gas flow 
lines, and always be in thermal equilibrium. Therefore 
the properties discussed for the carrier gas should be very 
similar for the species gas as well. If the species has inter- 
nal structure (i.e. electonic, vibrational, and rotational 
states that can be thermally populated) , then if there are 
sufficiently large inelastic, internal-state-changing colli- 
sion cross sections with the carrier gas then the internal 
temperatures can be thermalized as wellfl]. This is often 
the case, so supersonic beams are useful tools for creat- 
ing beams of translationally and internally cold molecular 
beams, albeit with a large forward velocity. 



E. Technical details of source construction 

The basic requirements for a buffer gas cell were out- 
lined in Section |II A[ but in this section we will provide 
some technical details of source construction. More de- 
tails may be found in existing literature [3l El [71 [9l |67] . 
Here we shall focus purely on the design and engineering 
aspects, so readers not concerned with these details are 
encouraged to skip this section. 



1. The cell 

Here we will describe the cell used by Hutzler et al.[7], 
however the features are similar for other ablation-loaded 



buffer gas cells [315115] (we will not consider the technical 
details of capillary loading [3 [67]). A schematic of the 
cell is shown in Figure [2] 

The cell is a cryogenically cooled, cylindrical copper 
cell with internal dimensions of 13 mm diameter and 75 
mm length. A 2 mm inner diameter copper tube enter- 
ing on one end of the cylinder flows buffer gas into the 
cell. A 150 mm length of the fill line tube is thermally 
anchored to the cell, ensuring that the buffer gas is cold 
before it flows into the cell volume. An open aperture 
(or nozzle) on the other end of the cell lets the buffer gas 
spray out as a beam. The aperture should be thin-walled 
(typically iO.5 mm) to prevent the species of interest from 
sticking to the sides of the aperture. The ablation target 
is located approximately 50 mm from the exit aperture. 
A pulsed Nd:YAG laser (Continuum Minilite II, 532 nm, 
17 mJ per 6 ns pulse), focused with a converging lens, is 
fired at the target. 

In addition to the main, cylindrical internal volume, 
two holes are drilled perpendicularly to the main vol- 
ume axis. One hole (13 mm diameter) has the ablation 
target on one side, and a window on the other; another 
hole (10 mm diameter) has a window on both sides for 
a spectroscopy laser. The windows are made of 3 mm 
thick borosilicate glass, and are sealed to the holes in the 
cell with indium. The window for the ablation laser is 
mounted on the end of a 38 mm long copper tube which 
attaches to the cell; this keeps the window far from the 
ablation target. If the window is too close to the abla- 
tion target, the window will become coated with ablation 
detritus, into which the ablation laser can begin deposit- 
ing energy. This both reduces the amount of energy de- 
posited into the ablation target, and cause the window 
to become damaged and unusable. 



2. Radiation shields 

The blackbody heat load on a 4 K surface from a 300 K 
environment is approximately 50 mW cm~^. A cryogenic 
cell can easily have surface area of >100 cm^, and would 
therefore have to experience a heat load of several watts. 
This would overwhelm most cryogenic refrigerators, so 
it is necessary to surround the cell with cold surfaces. 
Surrounding the cell by a blackbody shield cooled by a 
liquid nitrogen cryostat (77 K), or the first stage of a 
pulse tube cooler (^ 30 — 50 K) will reduce the blackbody 
heat load on the cell to < 1% of the 300 K blackbody 
heat load. This heat load will instead be deposited into 
the blackbody shield; however, a typical liquid nitrogen 
cryostat or pulse tube cooler can easily handle these heat 
loads at these temperatures. The blackbody heat load 
can also be further reduced by covering the blackbody 
shield with several layers of thin aluminum-coated mylar 
"superinsulation" . 

The cell is also sometimes additionally enclosed in an 
inner radiation shield, typically around 4-5 K. This can 
help reduce the heat load further, since it could be eas- 
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ier to superinsulate a radiation shield instead of the cell 
itself. Additionally, the 4 K shield could help to protect 
cryopumps, as discussed below. 



3. Cryopumps 

Enclosing the cell in radiation shields means that the 
gas conductance to the external vacuum chamber is very 
low. If no pumping occurred inside the radiation shields, 
the buffer gas would build up to the point where the 
residual pressure was too high to allow for a beam. In the 
case of neon buffer gas, this can be solved by surrounding 
the cell with a radiation shield kept at a temperature 
where the vapor pressure of neon is acceptably low, which 
is simple to achieve with standard cryogenic refrigerators 
or liquid helium (for example, the vapor pressure of neon 
is < 10~^ torr at 7 K). A neon cooled buffer gas beam 
inside a 4 K radiation shield can operate for extended 
periods of time; Hutzler et al.[7] operated such a beam 
continuously for over 24 hours without beam or vacuum 
degradation. 

For helium this is not practical; the radiation shield 
would have to be cooled to < 0.5 K for similar vapor pres- 
sures, which would require much more complex cryogenic 
systems and superfluid film management. Fortunately 
there is a solution, which is to use activated charcoal as 
a crvopump[55]. When cooled to < 10 K, activated char- 
coal becomes a cryopump for helium with up to several 
1 s~^ pumping speed per cm^ of charcoal, and can hold 
almost 1 STP liter of helium per gram (though these 
values are highly dependent on temperature, and other 
parameters [58]). Covering the inside of a 4 K radiation 
shield with ^ 100 cm^ of activated charcoal bonded to 
copper plates can result in a cryopump with a pump- 
ing speed of several hundred liters per second that can 
last for hours under typical gas loads [H [9]. As the char- 
coal fills up with helium the pumping speed begins to 
change, which can result in a degradation of beam prop- 
erties [21 [9]; therefore, it is necessary to warm up the 
charcoal (typically to ~ 20 K) and pump out the des- 
orbed helium with a mechanical pump. Experience shows 
that the beam properties are very sensitive to charcoal 
amount and placement [7j [9] , often in very non-intuitive 
ways. 



4- Beam collimation 

Similar to supersonic beams, buffer gas cooled beams 
typically have collimators to create differentially pumped 
regions and reduce the gas load on the part of the appa- 
ratus where the molecules are being probed. With buffer 
gas cooled beams this is especially advantageous, since 
collimators can be kept in the cryogenic region and there- 
fore take advantage of the very large pumping speeds af- 
forded by cryopumping. Hutzler et al.[7] found that the 
beam properties were largely insensitive to the placement 



of a collimator when using a neon buffer gas; however, as 
with the charcoal cryopumps, the placement of the colli- 
mator has often been seen to strongly influence the beam 
properties of the helium buffer gas cooled beam. Cor- 
rect placement requires careful consideration, and may 
involve multiple attempts. 



III. APPLICATIONS 

In this section, we will discuss some select applications 
where buffer gas beams can offer significant advantages: 
Laser cooling, precision measurements, collisional stud- 
ies, and trap loading. 

A. Laser cooling of buffer gas beams 

Laser cooling works by continuously scattering photons 
off of an atom or molecule in a manner that dissipates 
energv[48j. The number of photon scattering events to 
stop an atom or a molecule is given by A'' = mv/hk, 
where m is the mass of the species, v is its velocity, and 
k = 27r/A where A is the wavelength of the laser light. 
For typical atoms or molecules moving at room temper- 
ature thermal velocities, this number is about 10'''. The 
electronic ground state of an atom may contain multi- 
ple hyperfine states, and off-resonance excitation some- 
times brings the atoms to different hyperfine states in 
the excited state, which can decay to an off-resonant (or 
"dark") hyperfine state in the electronic ground state. If 
this happens, a "re-pump" laser, resonant with the other 
hyperfine state, can pump the atom back into the reso- 
nant state. 

In contrast, the electronic ground state of molecules 
contains multiple internal states, including vibrational 
and rotational, as well as hyperfine states. The num- 
ber of re-pump lasers needed to close the optical transi- 
tion grows accordingly. Despite these challenges, Di Rosa 
j93) pointed out that a set of molecules with strong sin- 
gle photon transitions and a highly diagonal vibrational 
transitions (Frank-Condon factors [85 ) may be feasible 
for direct laser cooling. The first criterion ensures a large 
photon scattering rate, meaning that a sufficient number 
of photons can be scattered in the limited interaction 
time (determined by the molecular beam velocity and 
geometrical constraints on the length of the interaction 
region). The second criterion implies a high probability 
of decaying back to the initial vibrational ground state 
[v" = 0). If p ^ 1 is the probability of decaying into a 
dark state after one optical transition, then, if there are 
no re-pump lasers, the total number of photons scattered 
for each molecule, N^cat, is given by (1 — p)^""* ^ e~^, 
yielding Nscat ^ 1/p- Di Rosa[n3] listed a few molecules 
that satisfied the above requirements, including CaH and 
NH. Most of these molecules can scatter ~ 100 photons 
without any vibrational re-pump lasers, and can scatter 
up to ^ 10** photons with one v" = 1 re-pump laser. 
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In contrast to vibrational leakage, where no strict se- 
lection rules apply, rotational transitions satisfy the total 
angular momentum change selection rule, A J — 0,±1. 
Stuhl et al. ^M] noticed that the rotational leakage can 
be mitigated by choosing a transition where the angu- 
lar momentum J" of the ground state is larger than the 
angular momentum J' of the lowest rotational state in 
the excited electronic state. For example, consider a 
molecule with a ground electronic state containing ro- 
tational levels J" = 1,2,3,.. ., and an excited electronic 
state with rotational levels J' = 0, 1, 2, . . . (for example, 
the X^Ai E^Uo transition of TiOl94 ). If a laser 
drives the transition \X, J" = 1) — > \E, J' — 0), then the 
selection rule \J" — J'\ < 1 means that the only radiative 
decay allowed is \E,J' = 0) -w \X,J" = 1), i.e. a ro- 
tationally closed transition. The molecule may still have 
vibrational leakage, but those channels can be closed with 
re-pump lasers in the usual way. 

Several of the diatomic molecules that satisfy the crite- 
ria for laser cooling as suggested by Di Rosa[n3] and Stuhl 
et al.[94] are radicals. Since the number of photons and 
the distance needed to stop molecules are proportional to 
W|| and respectively, where is the forward velocity of 
the molecules, molecular radical sources with a low for- 
ward velocity, as well as a high brightness, are desirable 
for direct laser cooling. As mentioned earlier (see Table 
|l]) , a buffer gas cooled beam can deliver higher brightness 
and lower forward velocity for these types of species. As 
an example, a boosted buffer gas beam has a typical for- 
ward velocity of ~ 150 m s~^, which reduces the number 
of photons needed to stop a CaH molecule by a factor 
of ~ 5 compared to a 1300 K oven source, as initially 
proposed by Di Rosa[93l . 

The first demonstration of optical cycling in molecules 
was achieved by using a buffer gas cooled beam of 
SrF^. The X^T,+^^ -> transition of SrF has 

relatively diagonal Franck-Condon factors [15] fv'=o~^v"- 
fo,o - 0.98, /o,i ~ 0.02, /o^2 « 4 X IQ-^, and /o,3+ < 
10~^. Theoretically, SrF could scatter Nscat = l//o,3+ > 
10^ photons with two vibrational re-pump lasers, assum- 
ing that other internal states are addressed as well. Rota- 
tional leaking was prevented by using an N" = 1 ^ N' ~ 
transition. Since SrF has a nuclear spin of / = 1/2 and 
the hyperfine splitting of |iV' = 0, J' = 1/2 > is small, ev- 
ery hyperfine state in the N" = 1 manifold would be pop- 
ulated during the transition and needs to be addressed. 
This was achieved using lasers with several sidebands 
(created using an electro-optic modulator) which were on 
resonance with the hyperfine states. One important fea- 
ture of a transition with J" > J' is that the ground state 
contains more Zeeman states than the excited state, lead- 
ing to dark Zeeman states [94] m,/" = ± J". To address 
these dark states, Shuman et al. used a small magnetic 
field which was at an angle relative to the polarization of 
the cooling lasers to constantly remix the Zeeman states. 

With the beam source used in those first laser cooling 
experiments, the SrF molecules were produced by laser 
ablation in a helium buffer gas cell at 4 K, resulting in a 



molecular beam with a total number of 10^ SrF molecules 
exiting the cell (per pulse) in the iV" = 1 state moving at 
V|| = 200 m s^^ [in]. The molecular beam was collimated 
transversely by a 6 mm diameter aperture before entering 
the laser interaction area, which consisted of the main 
cooling laser, one re-pump laser for w" — 1, and mag- 
netic field coils. When the magnetic field was turned on 
in addition to the main cooling laser, Shuman et al. ob- 
served that the induced fluorescence signals increased by 
a factor of ~ 3.5. Adding the vibrational re-pump laser 
further increased the fluorescence signals by another fac- 
tor of ~ 3.5, corresponding to a total average number 
of photons scattered per molecule of Ngcat 170. Since 
both lasers were applied perpendicularly to the molecu- 
lar beam from one transverse direction, the net momen- 
tum transferred from the photons to SrF was sufficient to 
create a position shift or deflection of the beam. The ob- 
served deflection indicated Ngcat = 140, consistent with 
the previous estimation based on the fluorescence mea- 
surement. 

In a subsequent experiment, Shuman et al. [50] demon- 
strated transverse cooling of the buffer gas cooled SrF 
beam. By incorporating an additional re-pump laser for 
the v" = 2 state and arranging the cooling lasers to in- 
tersect the molecular beam multiple times at nearly 90°, 
the molecules scattered Ngcat = 500 — 1000 photons in 
a 15 cm long region. When the cooling lasers were red 
detuned with a magnetic field of 5 G, the transverse tem- 
perature of the SrF beam was cooled from Tq = 50 mK 
(due to the geometry of the mechanical beam collima- 
tors) down to T ^ 5 mK. In this regime, the mixing 
field is large enough for rapidly remixing molecules in 
the Zeeman states, leading to Doppler cooling. Shuman 
et al. discovered that the molecular beam also obtained 
a lower temperature when the lasers were blue detuned 
with a relatively low field of 0.6 G. In this case, SrF in 
the bright Zeeman state constantly experienced a peri- 
odic potential provided by the cooling lasers, resulting in 
"Sisyphus" -like coohngdSHSU]. 



1. Direct loading of exotic atoms into magneto- optical traps 

Several exotic atoms have been successfully loaded 
into magneto-optical traps (MOTs[lSl[HS]), including Yb 
[HilSIl, Ag[nB], Er[M], Dy (TUD] and Tm [Ml. Because 
these atoms have much higher melting points than alkali 
atoms, effusive oven sources used to create these atomic 
beams operated at high temperatures, between 700 and 
1500 K. A Zeeman slower is typically needed to deceler- 
ate the fast-moving atoms (uj| ~ 350 — 580 m s^^) for 
efficiently loading into MOT, which has a typic al cap- 
ture velocity < 60 m s~^. As discussed in Section II C 2 



a near effusive buffer gas cooled beam can be produced 
by adding one slowing cell to the typical buffer gas cell. 
Although this type of beam has been found to have a 
typically lower extraction efficiency(~ 1%) than a hydro- 
dynamically extracted beam without a slowing cell, the 
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low forward velocity of the beam makes direct loading of 
exotic atomic beams into MOTs feasible. Additionally, 
the very high fluxes possible with ablation of certain met- 
als can lead to a superior beam system. 

The longitudinal and transverse temperatures of CaH 
molecules in a beam produced with a slowing cell were 
measured to be 5 K and ^ 2.5 K, respectively pD] , 
Since the cooling mechanism only relies on elastic col- 
lisions with the He buffer gas, we expect atomic beams 
produced from such a slowing cell to have properties simi- 
lar to the beam of CaH (except, of course, for any possible 
mass-related differences). For example, the forward ve- 
locity of a slow Yb beam with a longitudinal temperature 
of 5 K is ^ 30 m s~^, which is consistent with the forward 
velocity of a Yb beam demonstrated by Patterson [102] 
using a similar cell design. Employing a one-dimensional 
laser cooling model, the capture velocity of a MOT is 
given by Vc = (|A| +T)/k [97], where A is the detuning 
of the cooling laser and Y is the natural linewidth. For 
Yb, assuming | A| = F = 27r x 28 MHz, the capture veloc- 
ity is 22 m s~^, implying a large fraction of a buffer gas 
cooled Yb beam can be loaded without any longitudinal 
laser slowing (see Table IV I. 

CoUisional heating from the buffer gas streaming 
through the MOT region could be a concern. Typically, 
the buffer density inside the slowing cell is on the order 
of 10"'^'* cm^^. Although it is difficult to measure the He 
beam profile, a good approximation is to assume that 
the He beam attains the same temperatures as the CaH 
molecular beam. Therefore, the coUisional heating from 
the buffer gas is negligible compared to the cooling power 
from the scattered photons starting from a distance of 
10 cm away from the aperture, where the He density is 
< 10^^ cm^"^. Any potential loss due to the divergence 
of the beam can be mitigated by employing transverse 
laser cooling to coUimate the atomic beam between the 
cell aperture and MOT. Assuming the transverse capture 
velocity is given by v^^t = F/fc, we can estimate the load- 
ing efficiency of a buffer gas cooled beam into MOTs by 
considering the velocity distributions that are within Vc.t 
transversely and longitudinally. 

Table IIVI shows some exotic atoms that have been 
loaded into MOTs from effusive oven sources, and have 
also been buffer gas cooled. Of these species, only Yb 
has been realized as a buffer gas beam so far [5]. Given 
the generality of the buffer gas beam technique, we be- 
lieve that producing cold and slow exotic atomic beams 
with beam properties similar to the CaH beam [10] is 
very likely. We can therefore estimate the loading effi- 
ciency of these atoms into MOTs from buffer gas beams, 
fload, by assuming that these proposed atomic beams 
would have the same longitudinal and transverse tem- 
peratures of the CaH beam, 5 K and 2.5 K respectively. 
Based on the ablation yields inside the buffer gas cell, 
Nabiation, reported in the past [4] and a typical extrac- 
tion efficiency of the slowing cell of fext ^ 5 x 10^^, we 
can estimate potential loaded numbers per ablation pulse 
Nioad ^ Nabiation X fext X fload, showu as the last column 



Atom 


trapped, exp 


fload 


%] Nioad [pulse-1] 


Yb im 


10^ 


21 


2 X 10^° 


Ag [98] 


3 X 10*^ 


1.7 


3 X 10^ 


Dy [TOD] 


2.5 X 10** 


33 


3 X 10^ 


Er [99] 


10^ 


39 


4 X 10* 


Tm HOl] 


7 X 10"* 


0.3 


3 X to** 



TABLE IV. Some exotic atoms which have been loaded into a 
MOT from an effusive oven source, and also buffer gas cooled. 
Ntrapped,exp IS the Humber of atoms that were loaded into the 
MOT in the respective experiments, fioad is the projected 
loading efficiency from a buffer gas beam source, and Nioad 
is the projected number of atoms which could be trapped in 
the MOT with a single pulse from a buffer gas source. A 
description of these estimates may be found in the text. 



in Table IV The estimated numbers of atoms loaded per 
pulse are high compared to the reported trapped num- 
bers. By loading 1 — 10 pulses, the expected atom num- 
bers (or densities) accumulated in MOTs could approach 
a regime where the radiation pressure and collisions be- 
tween ground and excited state atoms would limit the 
atomic densitv |103j . Co- loading multiple species with- 
out using high temperature ovens and Zeeman slowers 
into a MOT is feasible by simply ablating several species 
simultaneously. 



B. Precision measurements 

Another application of molecular beams is the preci- 
sion determination of molecular properties. This is made 
possible by the lack of collisions in the beam region, 
and enhanced by the large volume of advanced molec- 
ular beam production, manipulation, and spectroscopy 
techniques [I] |2]. Precision spectroscopy with supersonic 
beams is a workhorse in chemistry and physics. As the 
reader can surely surmise by now, many precision mea- 
surements can benefit from the high flux and low for- 
ward velocity afforded by buffer gas cooled beams. In 
this section we will review some molecular beam preci- 
sion measurements, and discuss how they could benefit 
from buffer gas cooled beam technology. 



Electnc dipole moments 



In 1950, Purcell and Ramsev |104| were the flrst to se- 
riously consider the possibility of a fundamental parti- 
cle (i.e. an electron or nucleon) possessing a permanent 
electric dipole moment (EDM). Before this time, it was 
widely (though not universally [TOS]) believed that a fun- 
damental particle could not possess such a moment be- 
cause it would violate parity (P) and time-reversal (T) 
invariance in the following manner: Since a dipole mo- 
ment is a vector quantity it must be aligned with the spin 
of the particle, which is the only intrinsic vector associ- 
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ated with a fundamental particle. However, an electric 
dipole moment is a polar vector, while spin is an axial 
(or pseudo-) vector, so these two quantities must behave 
oppositely under inversion of coordinate systems or rever- 
sal of time, yet the previous argument requires them to 
uniquely determine each other. Although no permanent 
EDM of a fundamental particle has ever been measured 
(despite over 50 years of experimental searches 106i ), 
such a discovery would have a profound impact on our 
understanding of particle physics and cosmology. For an 
accessible discussion about EDMs of fundamental par- 
ticles and their implications, see the article by Fortson, 
Sandars, and Barr [107j : for more details, see the book by 
Khriplovich and Lamoreaux|28j . 

The essential idea of nearly every EDM search is 
to look for a small shift in precession frequency of a 
spin in an electric field as the electric field direction is 
reversed I28j. For a particle with spin S, magnetic dipole 
moment n and electric dipole moment d in a parallel 
magnetic field B and electric field £, the spin of the par- 
ticle (which is parallel to both the magnetic and electric 
dipole moment) will precess about the fields with angular 
frequency 



^ 'W\' 



(57) 



Consider an experiment where S and B are kept fixed 
but the direction of £ is reversed, and then the frequency 
uj± is measured for the two different orientations ±£. 
The value of the electric dipole moment d can then be 
extracted by comparing these two measured values. 



/i(w_i_ — w_) 
2£-S/\S\ 



(58) 



The determination of d then reduces to a precision mea- 
surement of the frequencies uj±. Precise frequency mea- 
surements arc best accomplished by the separated oscilla- 
tory field (or Ramsey) method[2], in which the spins are 
allowed to precess for a time (the "coherence time") t 
before the accumulated phase angle is measured, thereby 
determining ui. If such an experiment is repeated with 
a total of N spins, the uncertainty in the EDM is given 
bydg 



Sd- 



£WNI\S\ 



(59) 



From this equation, we can see that an EDM ex- 
periment should have large electric fields, long coher- 
ence times, and large count rates. Neutral objects such 
as neutrons or atoms can simply be placed in large 
electric fields; however; charged particles like electrons 
or protons would simply accelerate out of the electric 
field. Sandars[108j realized that not only could an elec- 
tric dipole moment of a charged object be inferred by 
studying a composite neutral object, such as an atom or 



molecule, but the electric dipole moment could actually 
be enhanced. For example, Sandars calculated ^09j that 
an electric dipole moment of the electron de could in- 
duce an electric dipole moment of an atom d^, with a 
relativistic "enhancement factor" dA/de that can be over 
100 for Cs atoms. The basic idea |110j is that the electron 
can sample the atomic electric field, which is typically of 
order 



£a 
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5 GV cm^ 



(60) 



where ao is the Bohr radius. This is considerably larger 
than a typical macroscopic lab field of < 100 kV cm~^. 
Atoms cannot take full advantage of this enhancement 
because they cannot be fully polarized in laboratory elec- 
tric fieldsflll]; however, polar molecules can be polarized 
in lab fields, and for this reason most current electron 
EDM experiments use polar molecules |111| . 

The molecules which have the highest enhancement 
factor have high-Z atoms with partially filled electron 
shells, and tend to be chemically reactive, have high melt- 
ing points, and are often free radicals [111) . For this rea- 
son, effusive beam sources of these molecules are typ- 
ically not feasible; the Boltzmann distribution of rota- 
tional states at temperatures where the species have ap- 
preciable vapor pressure put only a tiny fraction in any 
single rotational state (however, performing an electron 
EDM experiment with molecules in a hot vapor cell is a 
D0ssibilitv |112| ). A molecular beam for an EDM search 
should therefore have low temperature and, in light of 
Eq. (59), a low forward velocity to allow the molecules 



to interact with the fields for a long time, giving large 
T. For the types of molecules under consideration, buffer 
gas cooled beams [3 [H] have considerable advantages 
over the standard supersonic beam 111]. As discussed ear- 
lier, the brightness for refractory or chemically reactive 
species can be very high with buffer gas beams, and they 
have considerably slower forward velocities. In Table [ij 
we can see that a buffer gas cooled ThO beam has a 
slower forward velocity and over 100 times the bright- 
ness compared to a demonstrated, state-of-the-art YbF 
supersonic beam, both of which are currently being used 
for an electron EDM search. These reasons make buffer 
gas cooled beams an attractive option for electron EDM 
searches, with some existing supersonic experiments con- 
sidering switching to buffer gas cooling[70]. 

Searches for the electron EDM are currently underway 
in molecular beams of YbF[17], ThO[SS], and WC[m]. 
as well as several other atomic, molecular, and solid state 



systems |1 11 



2. Parity violation and anapole moments 

The exchange of weak neutral currents between the 
constituent particles of an atom or molecule can lead to 
observable parity violating effects in their spectra. Pre- 
cision examinations of these effects, known as atomic 
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parity violation (APV), allow the study of nuclear and 
electroweak physics using atomic systems [29]. An ex- 
perimentally feasible approach to using APV to probe 
electroweak physics was first proposed in 1974. |ll4] and 
since then APV has been observed in a number of 
atoms [29l I115H117| . Here we will briefly discuss how 
buffer gas beams may improve experimental studies of 
APV; detailed discussions of APV theory and experi- 
ments may be found elsewhere 1^ [1151 1118j . 

Weak neutral currents. The exchange of a neutral 
Zq boson between the nucleus and the electrons in an 
atom can give rise a parity-violating, electronic contact 
potential |114j . This potential can mix atomic states of 
opposite parity, such as a valence electron in ^1/2 and 
Pi/2 atomic orbitals. This mixing can either be mea- 
sured by searching for parity violating optical rotation 
in atomic emission/absorption |l 14|, I117j . or by looking 
for interference between the APV induced state mixing 
and the mixing provided by an external electric field (the 
Stark-interference method [119j ). APV has been mea- 
sured in a number of atoms |29l I116j , and comparison to 
theory gives tests of the standard model at low energies. 

Nuclear anapole moments. Another parity violating 
effect in an atom is the anapole moment 11181 1120j . 
which arises from an exchange of a Zq or W± boson be- 
tween nucleons in a nucleus. Like the parity-violating 
potential above, the interaction is a contact potential be- 
tween the electrons and nuclear anapole moment. The 
observable effect is very similar to that of the weak neu- 
tral currents discussed above, with the important ex- 
ception that the size of the effect will vary depending 
on the nuclear spin. An experimental signature of an 
anapole moment is then a difference in parity violat- 
ing amplitudes between different hyperfine levels. So 
far only a single anapole moment has been measured, 
in the ^■^•^Cs nucleus [1 2 IJ, though other experiments are 
underway |30j. 

Enhancement of parity violation in diatomic molecules. 
Parity-violating effects have thus far been observed only 
in atoms, though the search has been extended to di- 
atomic molecules for very compelling reasons. Molecules 
have states of opposite parity (rotational states, and in 
some cases A/r2-doublets[5S]) which are typically much 
smaller than the spacing between opposite parity elec- 
tronic states in atoms, so the effects of nuclear-spin- 
dependent parity-violating violation (such as anapole 
moments) can be greatly enhanced [T^IHUS]. This mix- 
ing can be increased even further by applying exter- 
nal magnetic fields to push rotational states to near- 
degeneracy. DeMille et al. |124| proposed (and are work- 
ing on) an experiment to use a supersonic beam of BaF 
free radicals in a electronic state to measure parity- 
violating amplitudes with high precision. 

Similar to the case with an EDM measurement, the 
shot-noise limited uncertainty (Eq. ( |59[ )) in the parity vi- 
olating amplitude scales as oc N~^^ T ^ , where N is the 
total number of interrogated molecules and T is the time 
that the molecules spend interacting with the electro- 



magnetic fields. Also similar to the EDM measurements 
is the choice of molecules, including ^30j BiO, BiS, HgF, 
LaO, LaS, LuO, LuS, PbF, and YbF: These are all free 
radicals with heavy nuclei and are therefore prime candi- 
dates for production in a buffer gas cooled beam. There 
are even some molecules which have already been dis- 
cussed, such as YbF and BaF (Table |l]), for which using 
a buffer gas beam can deliver over 100 times the bright- 
ness of a supersonic beam, with a slower forward velocity 
(allowing for the same time T with a shorter apparatus, 
easing technical requirements). It is for this reason that 
a BaF anapole moment experiment (at Yale |124| ') is pur- 
suing a buffer gas cooled beam source [T^. 

3. Time-variation of fundamental constants 

There has been much recent interest in the ques- 
tion of whether or not fundamental constants are truly 
"constant," or whether their values have changed over 
time. In particular, the dimensionless fine structure 
constant a = e^/47reo /ic and electron-proton mass ra- 
tio fi = m^/mp have attracted special attention due to 
the possibility of measuring their variation from mul- 
tiple independent sources. Detailed discussions of the 
theory and experiments discussed here may be found in 
reviews [551 [T^. 

Searches for time-variation of fundamental constants 
typically take one of two approaches. One is to look 
at data which may give information about a ot from 
a very long time ago, for example by examining astro- 
nomical spectra at large redshifts [551 1127j . but with lim- 
ited accuracy (though the search may be coupled with a 
laboratory precision spectroscopy experiment). Another 
technique is to use measurements with high precision but 
over a comparatively short timescale, for example by us- 
ing an atomic clock [128l I129j . The current limits on /i 
and a are 

^ = (1 ± 3) X 10-16 yr-i 

from the inversion spectrum of ammonia measured in 
sources with high redshift [127j , and 

- = (-1.6 ±2.3) X 10"^^ yr-i 
a 

from comparing optical clocks [T28] . 

As with EDMs and parity violation, molecules can have 
large enhancement factors which could allow for more 
sensitive measurements. In particular, the a and jj, de- 
pendence of fine structure splitting ujf cx a^, rotational 
splitting ljJj. (X ^ and vibrational splitting ijj^ oc /i^^^ are 
different [3S], and if two levels with different hyperfine, 
rotational, or vibrational character are nearly degener- 
ate there can be significant sensitivity enhancement [130] 
in changes to both and a. This enhancement can be 
several orders of magnitude (compared to atoms), and 
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there are several promising species including CuS, 
IrC, LaS, LaO, LuO, SiBr, YbF, Csz, and Si2- Several of 
these molecules may be difficult to produce through nor- 
mal beam methods, but are prime candidates for buffer 
gas beam production. 



C. Cold collisions and trap loading 

The study of cold molecular collisions may impact ar- 
eas of fundamental physics, astrophysics, and chemistry. 
For example, several cold neutral-neutral reactions are 
calculated to proceed with a large rate |131H133 \ and 
hence may play crucial roles in cold chemistry of inter- 
stellar clouds, where low mass stars may be formed [TM]. 
The number of direct measurements of reaction rates in- 
volving cold molecules is small compared to the number 
of theoretical calculations, but new measurements could 
be used in detailed astrophysical modeling, and can assist 
further theoretical calculations [iUl H^Tl 1132) . 

Collision studies at high temperatures may involve 
many internal states and scattering channels, which can 
complicate both theory and measurements; cooling the 
sample to a few K can reduce the involved states and 
channels, making the problem more tractable. In addi- 
tion, the potential energy provided by electromagnetic 
fields on magnetic or polar molecules (typically a few 
Tesla or kV cm~^, respectively) can be comparable to 
the thermal energy of the cold molecule, which makes 
the experimental control of collisions and chemical reac- 
tions an achievable goal |261 

Understanding cold molecular collisions is an impor- 
tant step for producing ultra-cold molecules via direct 
cooling methods. Cold molecules produced from direct 
methods (buffer gas cooling, supersonic expansion with 
Stark deceleration, etc.) so far possess a thermal en- 
ergy > fc^ X 1 mK. There are a few possible paths to 
approach ultra-cold temperatures, including evaporative 
cooling and sympathetic cooling of molecules. The for- 
mer (latter) relies on removing molecules (atoms) from 
the high energy tail in a thermal distribution, followed 
by subsequent elastic molecule-molecule (atom-molecule) 
collisions which rethermalize molecules to a lower tem- 
perature. Therefore, good molecule-molecule or atom- 
molecule collision properties are required for bringing 
cold molecules into ultra-cold temperatures. 

Buffer gas cooling can provide a large molecular sam- 
ple in a single quantum state with a kinetic energy on 
the order of ~ fc^x K. Several cold atom-molecule and 
molecule-molecule systems have been investigated using 
buffer gas cooling, including collisions of He atoms with 
2E and 3S molecules [Ml [135], N + NH[T36], and ND3 + 
OH[60 . 

The He -I- molecule and N + NH collision systems in- 
volved magnetic trapping of molecules and resulted in 
deeper understanding of atom-molecule interactions. For 
example, the He -I- and '^E molecule collision systems 
revealed the spin relaxation mechanisms of and '^E 



molecules in magnetic traps [1]. Cold N -I- NH collisions 
indicated not only that N might be a good candidate for 
sympathetic cooling of NH, but that a cold chemical re- 
action N -I- NH — >■ N2 -I- H may be suppressed in the 
magnetic field [136| . 



1. ND3 + OH collisions 

Trapping molecules provides an ideal environment for 
studying cold molecular collisions by providing a long in- 
teraction time. Using trapped molecules as the scattering 
target and monitoring the molecular trap loss allows ab- 
solute measurements of the collision cross sections |137j . 
Although several chemically interesting molecules have 
been cooled via direct methods and trapped, inelastic 
molecular collisions or reactions have not been observed 
in the trap because of a low molecular density and a low 
coUisional energy (~ 100 mK). Assuming a typical elastic 
cross section of 10^ ^^cm^, an elastic NH -I- NH collision 
occurs in hundreds of seconds with a molecular density 
of 10* cm"-^ at 0.5 K, which is much longer than the trap 
lifetime of ^1 s in the buffer gas loading experiment |135) . 

Absolute atom-molecule cross section measurements 
were obtained for the OH-I-D2 and OH-|-He systems by 
scattering supersonic D2 and He beams off of magneti- 
cally trapped OH molecules [137]. The lowest coUisional 
energies (84 K for He and 200 K for D2) achieved with 
a cooled supersonic nozzle were two order of magnitudes 
higher than the trapping potential of the molecules. In 
such a situation it is difficult to differentiate the contri- 
butions of elastic and inelastic collisions to the total cross 
section and control the collisions with an external field. 
Since buffer gas cooled beams of molecules can be cold, 
slow, and bright, they are an attractive choice for the 
study of molecular collisions. In addition, several buffer 
gas beams can provide continuous output (see Tabic |ll| , 
which could allow molecular interaction times compara- 
ble to the lifetime of the trapped species. 

Such an experiment was performed by Sawyer et al. 
[SU' by combining a buffer gas cooled beam of ND3 with 
trapped OH molecules. Collisions were studied on the 
~1 s timescale, and led to the first observation of cold 
(5 K) heteromolecular dipolar collisions. The continu- 
ous ND3 beam was guided to suppress the background 
buffer gas, giving a typical molecular density of 10* cm~^ 
and a mean velocity of 100 m s~^ before colliding with 
OH. A supersonic OH beam was Stark-decelerated and 
loaded to a magnetic trap with a density of ~ 10^ cm~^ 
at a temperature of 70 mK. To study the ND3 -I- OH 
collisions, the trap lifetime of OH in the presence of an 
ND3 beam was compared to the lifetime solely limited by 
collisions with the background gas, yielding a trap loss 
cross section of (Tg°y = 0.27 x 10~^^ cm^. At a colhsion 
energy of few K, only tens of partial waves participate 
in the collision, and hence external fields can affect the 
coUisional properties. After applying a large electric field 
to polarize both OH and ND3 , 0-'°^ increased by a factor 
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of 1.4. 

Sawyer et al. also performed both quantum scatter- 
ing and semi-classical calculations to compute the elastic 
and inelastic OH -I- ND3 cross sections. The calculated 
cross sections were orders of magnitude larger than typi- 
cal gas-kinetic values, possibly resulting from the dipolar 
interaction. To calculate the cross sections in the field, 
Sawyer et al. included a long range and anisotropic dipo- 
lar interaction term C^/r^, in addition to an isotropic 
van der Waals interaction for the zero field case. They 
discovered that the elastic cross section increased by a 
factor of 2.7 in the presence of a field, compared to a 
relatively constant inelastic cross section. After consid- 
ering the probability of an elastic scattering leaving OH 
un-trapped. Sawyer et al. found the value of (7^°^^ and 
its increase in the field were consistent the theoretical 
model, indicating the elastic OH + ND3 collisions not 
only dominate the trap loss but can be controlled by the 
external electric field. 

Compared to other beam technologies, buffer gas 
cooled beams are well suited for these types of colli- 
sion experiments. Stark-deceleration of a supersonic ND3 
source can yield a beam with a forward velocity of ~ 90 m 
s~^ and a peak density of ~ 10* cm~'^, but with pulsed 
operation, < 1 ms pulse duration, and repetition fre- 
quency of few Hz[TB]. The fractional loss of trapped 
molecules due to colliding with another beam is given 
by 



.a'^'^VrelAt, (61) 



where nbeam is the density of the beam and Vrei is the 
mean relative velocity. Using a continuous buffer gas 
cooled ND3 beam increases the average interaction time 
by a factor of '^lOO compared to a Stark-decelerated ND3 
beam, leading to a higher sensitivity in observing cold 
collisions. 



2. Direct trap loading of molecules 

Cold and slow beams of CaH, with a kinetic energy 
of few fcfexK, can be created by using a slowing cell (see 
Section II C 2 1 , which should allow direct loading of the 



molecules into a magnetic trap. A magnetic lens or guide 
can be used to ensure that the buffer gas atoms do not 
interfere with the trapped molecules, similar to the ex- 
periment of Sawyer et al.[Sn] discussed earlier. Here we 
will explore one possible loading scheme, which combines 
magnetic slowing and optical pumping to achieve trap 
loading of slow CaH molecules. By using a buffer gas 
cooled beam, the proposed method benefits from large 
molecular fluxes, chemical versatility, and high vacuum 
in the trapping area. A similar scheme had been em- 
ployed to continuously load an atomic Cr beam into an 
optical dipole trap[138L 

Figure [I6] illustrates the scheme to continuously load 
CaH molecules into a quadrupole trap. When a low-field 
seeking molecule (m" = 1/2) enters the trap, it travels up 



the potential hill and slows down. Near the saddle point, 
a laser can optically pump the molecules to the high- 
field seeking state (m" — —1/2). Molecules continue to 
be decelerated while approaching the center of the trap, 
when another laser pumps the molecule into the trapped 
m'j —1/2 state. In this loading scheme, molecules with 
a kinetic energy large enough to overcome two potential 
hills can reach the trap center. Molecules with a kinetic 
energy lower than the trap depth after deceleration will 
remain trapped. For a magnetic quadrupole trap with 
a depth of 4 T, the capture energy for CaH is 5-7 K, 
suggesting that a large fraction of the slow CaH beam 
can be captured. 

To reduce the background buffer gas density in trap 
area, a magnetic lens can be used to collimate the 
molecular beam while the buffer gas is pumped away 
with charcoal cryopumps 



II El. We per- 



(see Section 

formed a semi-classical Monte Carlo simulation to cal- 
culate the efficiency of loading the CaH beam into the 4 
T deep quadrupole trap, indicating a typical efficiency of 
^ 5 X 10"'^, or about 3 x 10^ trapped molecules per pulse. 
The background gas density in the trapping region can 
be kept low such that molecules can be accumulated in 
the trap via continuously loading, leading to a total num- 
ber of 10'^ - 10* over time. This compares favorably to 
the typical values of 10"'"^^ molecules loaded into a trap 
from a Stark-decelerated supersonic source [T71 1139I - [T^ . 

We would like to point out that this scheme is gen- 
eral and applicable to other magnetic species including 
atoms. More importantly, it only requires scattering a 
few photons to achieve loading and therefore does not 
require a highly closed cycling transition. This opens 
up the possibility of trapping a large class of magnetic 
molecules. 



n,„ (j'=1/2, m=+1/2) 




m,=-Fl/2 



m =-1/2 



i:+(j"=i/2) 



FIG. 16. Loading scheme of a slow buffer gas beam into a 
magnetic trap. Pink and black curves indicate the poten- 
tial energy of the low-field and high-field seeking states in a 
magnetic trap, respectively. Red and blue arrows show the 
optical pumping processes to achieve irreversible loading due 
to spontaneous emissions. 
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